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NOTATION 


a  =  minimum  value  for  the  beta  distribution 
Ao  =  original  area 

Ap  =  area  of  the  piston  in  the  triaxial  extension  test  apparatus 
b  =  maximum  value  for  the  beta  distribution 
C  =  seismic  compression  wave  velocity 
cov  =  covariance 

dV/Vo(cyl)  =  volumetric  strain  calculated  assuming  deformation  as  a 
right  circular  cylinder 

E  =  expected  value 

F  =  ratio  of  the  between  sample  variance  and  the  within  sample 
variance 

F  =  area  under  the  right-hand  tail  of  the  F  distribution 
a 

Fq  q5  =  area  under  the  right-hand  tail  of  the  F  distribution  for 
*  a  =  0.05 

Gs  =>  specific  gravity 

HQ  =  original  height 

HC  =  hydrostatic  compression 

L  =  load 

LC  =  lower  chamber  pressure  in  the  triaxial  extension  test 
apparatus 

M  =  constrained  modulus 
Ml  =  loading  constrained  modulus 
My  =  unloading  constrained  modulus 
MNS  =  mean  normal  stress 

n  =  window  size  in  the  smoothing  technique  and  total  number  of 
curves  in  a  set  of  data 

PSD  =  principal  stress  difference 

r  =  sample  correlation  coefficient 

S  =  degree  of  saturation 

S-| ,  S2  =  secant  slopes 

Sxx  =  linear  regression  coefficient  for  the  x-axis  values 

Syy  =  linear  regression  coefficient  for  the  y-axis  values 

Sxy  =  linear  regression  coefficient  for  the  x-  and  y-axis  values 

TXC  =  triaxial  compression 

TXE  =  triaxial  extension 
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UC  =  upper  chamber  pressure  in  the  triaxial  extension  test 
apparatus 


B 


UX  = 
V  = 

Va  = 
x  = 

= 

ys(Xi)  = 
a  = 
B  = 

STRAIN  = 

y  = 

Yd  = 
AH  = 

Ax  = 

Ay  - 


'IN  SITU 
£s  = 


Gz 

VI  = 

\  = 
JUN  = 


UN 


ac  = 


z 

<p  = 
0)  = 
R  = 


uniaxial  strain 
variance 

air  voids  content 
mean  value 
measured  data 
smoothed  data 

beta  distribution  parameter 
beta  distribution  parameter 

ratio  of  in  situ  axial  strain  to  seismic  axial  strain 
wet  density 
dry  density 

change  in  height  of  a  test  specimen 
increment  size  along  the  x-axis 
increment  size  along  the  y-axis 
in  situ  axial  strain 
seismic  axial  strain 
radial  strain 
true  radial  strain 

axial  strain 

number  of  increments  in  the  covariance  analysis 
loading  Poisson's  ratio 
unloading  Poisson's  ratio 
initial  unloading  Poisson’s  ratio 

correlation  coefficient 
standard  deviation 
confining  pressure 
radial  stress 
axial  stress 

angle  of  internal  friction 
water  content 
incremental  slope 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


_ Multiply _ _ 

degrees  (angle) 

feet 

inches 

kips  (force) 

kips  (force)  per  square  inch 

megatons  (nuclear  equivalent 
of  TNT) 

pounds  (force)  per  square  inch 
pounds  (mass) 

pounds  (mass)  per  cubic  foot 


By 

0. 01745329 
0.3048 
25.4 

4.448222 

6.894757 

4.184 


_ To  Obtain 

radians 

metres 

millimetres 

kilonewtons 

megapascals 

petaj oules 


0.006894757  megapascals 

0.4535924  kilograms 

16.01846  kilograms  per  cubic 

metre 
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A  TECHNIQUE  FOR  DEVELOPING  PROBABILISTIC 


PROPERTIES  OF  EARTH  MATERIALS 


CHAPTER  1 
INTRODUCTION 


1.1  BACKGROUND 

The  vulnerability  of  buried  structures,  such  as  missile  silos  and 
shelters,  to  nuclear  blasts  is  assessed  by  conducting  scaled  test  events 
using  conventional  explosives.  The  stress,  velocity,  and  deflection  of 
the  structure  and  surrounding  soil  are  measured.  The  response  of  the 
buried  structure  is  signif icantly  affected  by  the  response  of  the  sur¬ 
rounding  earth  media.  The  stress  being  transmitted  to  the  structure  by 
the  soil,  the  friction  at  the  soil-3tructure  interface,  the  motion  of 
the  soil  with  respect  to  the  structure,  and  the  cratering  characteris¬ 
tics  of  the  soil  are  all  important  in  the  survivability  of  buried 
structures.  To  effectively  predict  or  evaluate  the  results  from  a  test 
event,  the  properties  of  the  earth  materials  must  be  determined.  This 
requires  extensive  soil  exploration  programs  involving  (a)  drilling  and 
sampling  to  obtain  disturbed  and  undisturbed  samples  for  identification 
and  laboratory  testing,  (b)  seismic  surveys  to  determine  the  seismic 
velocity,  (c)  mass-volume,  index,  strength,  and  compressibility  tests  on 
disturbed  and  undisturbed  samples,  and  (d)  in  situ  compressibility  tests 
using  explosives  as  the  loading  source. 

The  ground  shock  calculation  techniques  currently  used  by  the 
analysts  of  explosive  test  events  require  deterministic  input  parameters 
and  provide  deterministic  results,  i.e.,  single-valued  quantities  or 
functions  are  input  and  output.  In  actuality,  the  properties  of  earth 
materials  and  the  characteristics  of  the  explosive  blasts  are  dispersed 
random  variables.  Since  these  input  variables  arc  random,  the  calcu¬ 
lated  stresses  and  ground  motions  are  also  random  variables.  Therefore, 
ground  shock  problems  should  be  analyzed  probabilistically.  By  treating 
a  ground  shock  calculation  probabilistically,  the  effect  of  the  uncer¬ 
tainties  in  the  material  properties  and  characteristics  of  the  explosive 


blast  on  the  dispersion  of  the  output  variables  can  be  evaluated. 
Ultimately,  the  probability  of  failure  for  a  buried  structure  can  be 
assessed  if  the  uncertainty  in  the  input  is  known.  To  this  end,  a 
three-phase  project  is  underway  at  the  US  Army  Engineer  Waterways  Exper¬ 
iment  Station  to  (a)  develop  a  procedure  for  providing  the  probabilistic 
properties  of  earth  materials  needed  as  input  for  a  probabilistic  ground 
shock  calculation,  (b)  develop  a  probabilistic  cap-type  elastic-plastic 
constitutive  model  that  describes  the  probabilistic  properties  of  the 
earth  materials  and  does  not  require  much  more  effort  to  fit  than  the 
associated  deterministic  model,  and  (c)  develop  a  probabilistic  ground 
shock  calculation  computer  code  and  calculation  procedure.  This  report 
documents  the  development  of  phase  (a). 

1.2  LITERATURE  REVIEW 

It  has  long  been  recognized  that  probabilistic  concepts  and  methods 
should  be  used  in  geotechnical  engineering  to  aid  the  development  of 
practical  procedures  dealing  with  risk  and  uncertainty.  Efforts  have 
been  undertaken  to  analyze  probabilistically  subgrade  soil  strengths 
(Reference  1),  settlement  (Reference  2),  and  slope  stability  (Refer¬ 
ences  3  and  H).  In  developing  basing  concepts  for  ballistic  missiles, 
it  was  recognized  that  probabilistic  analysis  of  ground  shock  tests  was 
a  necessity  (Reference  5).  The  original  basing  concepts  for  the  MX 
missile  essentially  consisted  of  probabilistic  "shell  games"  in  which  a 
relatively  small  number  of  missiles  were  concealed  within  a  much  larger 
number  of  hardened  shelters.  If  the  nuclear  ground  shock  calculations 
used  in  shelter  vulnerability  analyses  were  to  be  consistent  with  this 
scenario,  it  was  necessary  that  the  soil  profiles  and  properties  used  as 
input  be  probabilistic  quantities  (Reference  6). 

Chou  used  a  concept  of  spatial  average  to  probabilistically  model 
subgrade  strength  for  pavement  design  (Reference  1).  He  pointed  out 
that  a  probability-based  design  procedure  would  allow  for  (a)  design  to 
a  level  of  reliability  depending  upon  the  importance  and  function  of  the 
pavement,  (b)  quantifying  the  design  risk,  and  (c)  optimizing  design 
results  by  controlling  and  minimizing  expenditures  and  early  failures. 
The  probabilistic  subgrado  strengths  used  in  the  procedure  were  based  on 
the  concept  of  spatial  averages  outlined  by  Vanmarcke  in  Reference  7. 
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The  application  of  spatial  averages  has  the  advantage  of  reducing  the 
variance  of  the  variables,  thereby  increasing  the  reliability  of  pre¬ 
dicted  pavement  performance. 

Diaz-Padilla  and  Vanmarcke  developed  a  probabilistic  soil-structure 
interaction  model  that  yields  first-order  probabilistic  information 
about  the  differential  settlements  of  a  foundation  in  terms  of  the 
probabilistic  input,  i.e.,  loads  and  soil  properties  (Reference  2).  The 
permanent  loads  vary  due  to  discrepancies  between  specified  and  actual 
member  dimensions  and  fluctuations  in  volumetric  weights.  The  dead 
loads  vary  since  the  final  location,  composition,  and  weight  of  non- 
structural  items  is  not  known  in  the  design  stage.  Live  loads  are 
random  since  they  are  impossible  to  predict  with  certainty.  The  mechan¬ 
ical  properties  of  an  assumed  nominally  homogeneous  soil  stratum  exnibit 
a  considerable  amount  of  variability  from  point  to  point.  The  random¬ 
ness  is  due  to  the  natural  variability  of  the  material  (the  mineral 
composition,  stratum  formation  process,  etc.)  and  uncertainty  introduced 
when  measuring  the  soil  properties  (sample  disturbance,  human  errors, 
and  test  imperfections).  The  means  and  variances  of  the  differential 
settlements  obtained  from  the  application  of  the  probabilistic  model  can 
be  used  to  estimate  the  probability  of  exceeding  specified  levels. 

Gilbert  developed  a  probabilistic  approach  for  analyzing  embankment 
stability  on  soft  saturated  clay  at  the  end  of  construction 
(Reference  3).  The  model  was  proposed  as  a  supplement  to  the  design 
approach  based  on  <j>  =  0  analyses  using  the  undrained  strength  of  the 
foundation  clay.  Two  case  studies  analyzed  with  the  proposed  model 
indicated  that  bias  inherent  in  field  vane  and  unconfined  compression 
testing  can  be  the  major  source  of  uncertainty  in  design. 

Vanmarcke  presented  a  probabilistic  approach  to  the  classical  limit 
equilibrium  slope  stability  problem  (Reference  4).  He  noted  that  risk 
analysis  should  be  U3ed  to  supplement  the  conventional  procedures  for 
determining  the  factor  of  safety  against  shear  failure  of  earth  slopes 
and  embankments.  Risk  analysis  permits  rational  consideration  of  fac¬ 
tors  that  significantly  influence  the  probability  of  failure,  such  as 
the  inherent  variability  of  soil  strength,  the  manner  in  which  strength 
values  are  selected,  and  the  amount  and  quality  of  soil  exploration  and 
testing. 


3 


Since  the  soil  parameters  have  a  significant  influence  on  stability 
calculations,  they  are  worth  quantifying.  Vanmarcke  noted  in  Refer¬ 
ence  7  that  the  variability  of  soil  properties  should  be  considered  in 
the  design  and  evaluation  of  soil  exploration  programs.  Three  major 
sources  of  uncertainty  in  modeling  a  soil  profile  are:  (a)  the  spatial 
in  situ  variability  of  the  soil,  (b)  statistical  uncertainty  due  to  a 
limited  number  of  soil  samples,  and  (c)  measurement  errors  due  to  sample 
disturbance,  test  imperfections,  and  human  factors.  The  first  source  of 
uncertainty  can  be  examined  by  considering  the  spatial  variation  of  the 
soil  properties.  The  effect  of  the  second  source  of  uncertainty  can  be 
reduced  by  conducting  many  tests.  The  effects  of  the  third  source  can 
be  reduced  by  careful  monitoring  to  obtain  high  quality  test  data. 

A  significant  use  of  probabilistic  ground  shock  analyses  is  to 
evaluate  the  variability  in  the  input  parameters  (blast  loading  and  soil 
properties)  in  terms  of  their  effects  on  the  dispersion  of  the  output 
variables  (stress  and  motion).  Using  the  method  of  partial  derivatives 
(Reference  8),  Rohani  (Reference  9)  developed  a  probabilistic  solution 
for  one-dimensional  (ID)  plane  wave  propagation  in  homogeneous  bilinear 
hysteretic  materials.  The  solution  was  incorporated  into  a  computer 
code  by  Rohani  and  Cargile  (Reference  10).  The  code  can  be  used  to 
quantitatively  rank  the  relative  effects  of  input  variabilities  on  the 
dispersion  of  the  output  quantities.  The  results  from  ID  probabilistic 
calculations  indicated  that,  at  early  times,  uncertainty  in  the  soil 
compressibility  had  the  greatest  effect  on  the  output  uncertainty.  The 
influence  of  soil  compressibility  uncertainty  decreased  with  increasing 
time  while  the  contribution  due  to  airblast  impulse  uncertainty  gradu¬ 
ally  increased. 

A  method  for  objectively  comparing  measurements  of  waveforms  from 
an  explosive  test  event  with  those  obtained  from  probabilistic  predic¬ 
tion  calculations  was  developed  by  Baladi  and  Barnes  (Reference  11).  By 
calculating  the  magnitude  and  phase-and-frequency  errors  between  the 
measured  and  expected  waveforms,  an  objective  judgment  of  the  degree  of 
agreement  or  disagreement  can  be  made. 

A  "simple"  method  was  used  to  develop  the  probabilistic  mechanical 
responses  presented  by  Vanmarcke,  Jackson,  and  Akers  (Reference  12)  and 
Jackson  (Reference  13).  The  method  involved  calculating  the  mean  and 
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standard  deviation  values  for  strain  at  given  values  of  stress  for  a  set 
of  uniaxial  strain  compressibility  tests. 

1.3  PURPOSE 

The  purpose  of  this  investigation  was  to  document  an  analysis  tech¬ 
nique  for  developing  probabilistic  properties  of  earth  materials.  The 
technique  involves  (a)  determining  the  subsurface  profile  for  a  given 
test  site  based  on  mass-volume  relations,  index  parameters,  seismic 
velocities,  material  type,  etc.,  and  (b)  developing  probabilistic  values 
or  relations  for  the  mass-volume  relations,  index  parameters,  seismic 
velocities,  and  mechanical  properties  (i.e.,  strength  and  compressibil¬ 
ity  relations)  for  each  layer  of  the  profile.  The  probabilistic  proper¬ 
ties  are  based  on  laboratory  test  results,  field  experiments,  and  engi¬ 
neering  knowledge.  Generally,  the  data  obtained  from  mechanical  prop¬ 
erty  tests  do  not  extend  to  the  stress  levels  required  for  ground  shock 
calculations.  A  procedure  for  extrapolating  the  variance  of  the  lower 
stress  portion  of  the  recommended  probabilistic  mechanical  properties  to 
higher  stresses  is  also  presented. 

1 .  SCOPE 

The  procedures  for  (a)  developing  the  subsurface  profile  of  a  given 
test  site  and  (b)  determining  the  probabilistic  values  and  relations  for 
the  mass-volume  relations,  index  parameters,  seismic  velocities,  and 
mechanical  properties  for  each  layer  of  the  profile  are  described  in 
Chapter  2,  A  technique  for  removing  the  electrical  noise  that  is  often 
embedded  in  the  measured  data  from  mechanical  property  tests  and  Rosen- 
blueth's  point  estimate  procedure  are  also  explained  in  Chapter  2.  The 
application  of  the  analysis  technique  to  data  obtained  from  two  adjacent 
areas  located  near  Yuma,  Arizona  is  presented  in  Chapter  3.  The  recom¬ 
mended  probabilistic  mechanical  properties  are  presented  in  Chapter  iJ. 

A  method  for  extrapolating  the  variance  of  the  lower  stress  mechanical 
properties  to  higher  stresses  is  also  presented  in  Chapter  4.  A  sum¬ 
mary,  conclusions,  and  recommendations  for  further  study  are  presented 
in  Chapter  5. 
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CHAPTER  2 


DESCRIPTION  OF  THE  ANALYSIS  TECHNIQUE 

2.1  DETERMINATION  OF  THE  SUBSURFACE  PROFILE 

The  subsurface  profile  of  an  area  is  based  on  changes  in  material 
type,  mass-volume  relations,  index  values,  seismic  velocity,  etc.  This 
data  is  plotted  versus  depth  to  determine  trends.  The  mean  and  standard 

. — ,  A 

deviation  values  (  x  and  o  ,  respectively)  are  calculated  for  a  zone 
of  data  which  appears  to  be  different  from  the  data  above  and  below 
it.  These  statistics  are  used  in  a  one-way  analysis  of  variance  (Refer¬ 
ence  14)  to  test  the  null  hypothesis  that  the  apparent  difference  in  the 
mean  values  for  adjacent  zones  is  not  statistically  significant. 

In  the  one-way  analysis  of  variance,  the  null  hypothesis  is  true  if 
the  ratio  of  the  between  sample  variance  and  the  within  sample  variance 

(  F  )  is  greater  than  F  ,  where 

a 


k 

SS(Tr)  iJ1  (m.  -  1) 
SSE  (k  -  1) 


and, 


SS(Tr)  = 


k 

I 

i-1 


m. 

l 


<*i> 


SSE  =  Z  (9  2  (m.  -  1)) 
i=1  1 


(2.1) 


(2.2) 


(2.3) 


In  equations  2. 1-2. 3,  k  is  the  number  of  samples  and  m  is  the  number 

of  data  in  a  given  sample.  The  value  for  F  is  obtained  from  the 

a 

k 

distribution  of  F  using  k-1  and  z  (m.  -  1)  degrees  of  freedom. 

i=1  1 

If  the  test  indicates  that  the  means  are  not  signil icantly  different, 

i.e.,  F  is  less  than  F  ,  the  zone  interfaces  are  redrawn  and  new 

a 

mean  and  standard  deviation  values  are  calculated.  This  process  is 
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repeated  until  the  one-way  analysis  of  variance  test  of  the  null  hypo¬ 
thesis  for  adjacent  zones  is  false,  i.e.,  F  is  greater  than  F^  . 

2.2  ANALYSIS  OF  MECHANICAL  PROPERTY  TEST  DATA 

Mechanical  property  tests  (i.e.,  strength  and  compressibility 
tests)  are  conducted  on  undisturbed  samples  from  each  layer  of  the 
subsurface  profile.  The  mechanical  property  test  results  for  samples 
from  a  given  layer  are  grouped  together,  electrical  noise  is  removed 
from  the  data,  and  a  covariance  analysis  is  conducted  on  the  data  set. 

2.2.1  Data  Smoothing 

Artificial  noise  created  by  the  electronics  of  the  measurement 
system  used  during  a  mechanical  property  test  is  embedded  in  the  test 
data.  For  some  tests,  this  noise  can  be  severe  and  almost  mask  the 
measured  response  of  the  specimen.  Therefore,  a  technique  to  smooth  the 
measured  data  without  changing  the  intrinsic  response  is  needed.  Such  a 
procedure  was  developed  by  Baladi  and  Barnes  (Reference  11)  and  was 
based  on  the  concept  of  a  marching  mean  square.  If  the  measured  value 
of  the  ith  data  point  is  expressed  as  y  (x^)  ,  the  corresponding 
smoothed  response  Ys(xi)  can  be  expressed  as 


ys(V  ■ 


IT 

n-1 


<V 


(2.4) 


where  n-1  is  the  window  over  which  the  marching  mean  square  is  taken 

(i.e.,  is  the  number  of  data  points  to  the  left  and  to  the  right 

of  the  ith  data  point).  The  value  of  n  must  be  an  odd  integer  equal 
to  or  greater  than  3-  The  degree  of  smoothing  increases  as  n 
increases.  A  small  value  for  n  should  be  used  initially.  The  data 
can  be  smoothed  repeatedly  with  the  same  window  size  until  a  satisfac¬ 
tory  response  is  produced  or  until  further  repetitions  produce  no  change 
in  the  "smoothed"  data.  If  additional  smoothing  is  required,  a  larger 
window  size  must  be  used.  The  smoothed  data  should  follow  the  general 
trends  of  the  original  data.  The  results  from  the  smoothing  process 


7 


should  be  compared  with  the  original  data  to  ensure  that  the  general 
trends  were  not  altered  by  excessive  smoothing. 

2.2.2  Covariance  Analysis  of  Smoothed  Data 

A  generic  procedure  for  statistically  analyzing  nonlinear  data  was 
developed  by  Baladi  and  Rohani  (Reference  15).  The  procedure  is  out¬ 
lined  in  this  section.  The  objective  of  the  analysis  is  to  determine 
the  mean  curve  with  its  one  standard  deviation  bound  for  a  set  of  curves 
by  relating  the  random  variables  y  and  x  (Figure  2.1).  This  can  be 
accomplished  by  applying  standard  statistical  procedures  to  the  slope  of 
the  random  curves.  The  following  steps  should  be  taken  to  conduct  the 
statistical  analysis: 

(1)  For  a  given  set  of  n  curves,  divide  the  x-axis  into  ji 
number  of  equal  increments  Ax  .  (This  procedure  may  also  be  applied  by 
dividing  the  y-axis  into  p  number  of  equal  increments  Ay  .) 

(2)  For  the  ith  increment,  determine  the  slope  of  the  jth 
curve  denoted  by  52^ 


Lyi 

fii j  =  Ax.  *  ^  ~  1 > 2» • • *n  (2.5) 

(3)  Determine  the  expected  value  and  the  standard  deviation  of  the 
slope  at  the  ith  increment  for  all  the  curves  according  to  the  follow¬ 
ing  expressions 


«1  =  E  (52.) 


i  2 

n  j,i 


52.  . 
1J 


(2.6) 


o  (52.) 


n 

l 

j=1 


(2.7) 


(4)  Next,  compute  the  mean  and  the  standard  deviation  of  y  .  To 

accomplish  this,  the  covariance  and  the  correlation  coefficient  matrices 

of  the  slopes  cov(52,  ,  52  )  and  0  ,  respectively,  are  first  calculated 

k  m 

from  the  following  relations: 
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cov  (n  , q  )  =  e  C(q.  -  n,,)(Qm  -  o  )3 

km  k  k  m  m 


(2.8) 


1 

n-1 


n 

E 

j=1 


(Skj" 


V 


(n  -  a > 

mj  m 


COV  ( Q  ,  fi  ) 
k  m 


km 


e  [(fl  -  a.)2]  e  [(a  -  a  )”] 

!/  v  ni  m 


m 


(2.9) 


In  which 


s  [(!,k  -  vS]  ■  ^  j,  %-v2  (2-io> 

where  k  =  1,2,..,!,..  ji  and  m  =  1,2,..,i,..  p  . 

Finally,  the  mean  value  and  the  standard  deviation  of  y  at  the 
i^h  increment  become 


_  i 

y.  =  E  E  [a.]  Ax. 
1  j=1  J  J 


(2.11) 
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This  procedure  is  well  suited  for  providing  probabilistic  mechani¬ 
cal  properties  because  the  slopes  of  the  plotted  responses  are  used  in 
the  analysis.  It  is  the  slopes  of  plotted  mechanical  responses  that  are 
fit  to  constitutive  models. 

The  covariance  technique  can  be  applied  in  two  ways.  First,  the 
technique  can  be  applied  to  smoothed  data  that  is  plotted  in  a  form  that 
is  conducive  to  constitutive  modeling.  Or,  the  technique  can  be  applied 
directly  to  the  smoothed  data,  and  the  results  of  the  analysis  converted 
into  a  form  that  is  conducive  to  constitutive  modeling.  A  procedure 
aeveloped  by  Rosenblueth  can  be  used  for  the  conversion.  These  two 
approaches  will  be  applied  and  compared  in  Chapter  3« 
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2.2.3  Rosenblueth’s  Point  Estimate  Procedure 


Rosenblueth  ceveloped  a  procedure  for  determining  the  moments  of  a 
dependent  variable  in  terms  of  functions  of  the  moments  of  its  indepen¬ 
dent  variables  (References  16  and  17).  The  usual  method  of  obtaining 
approximate  formulas  from  the  Taylor  expansion  of  a  function  about  the 
expectations  of  the  random  variables  (Reference  8)  is  bound  by  excessive 
restrictions  on  the  function,  such  as  existence  and  continuity  of  the 
first  few  derivatives.  The  Rosenblueth  procedure  avoids  such  restric¬ 
tions  by  using  point  estimates  of  the  function,  which  leads  to  expres¬ 
sions  similar  to  finite  differences. 

If  Y  is  a  function  of  two  random  variables  X1  and  X2  » 

Y  =  Y  (x1#  X2) 

and  the  probability  distribution  functions  of  X-|  and  X2  are  symmet¬ 
rical,  the  expected  value  of  Y  ,  E(Y)  ,  and  variance  of  Y  ,  V(Y)  , 
can  De  estimated  by  the  following  expressions: 


E(Y)  =  Y++  +  Y+_  (— Y“)  ♦  Y'+  (  — +  Y"  j  (2.13) 


V  (Y)  =  E  (Y2)  -  E2(Y) 


(2.111) 


where , 


E(Y2) 


++ . 


^12\ 

__  2  1‘LP1 


(Y  )  ( — f.— I  +  (Y  ) 
J1  2 


_+  2/1_P 


12 


(2.15) 


+  (Y  ) 


4 


and 


Y±±  =  Y  (X1  ±  5  X2  1  °X2} 


The  terms  x^  snd  x  are  the  expected  values  of  the  variables  X^ 

and  X2  ,  respectively;  likewise,  o  and  o  are  the  standard 

X1  2 
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is  the  correlation 


deviations  of  those  random  variables.  The  term 
coefficient  of  X-|  and  ^  • 

From  the  above  equations,  the  expected  value  and  the  variance  of 
Y  can  be  calculated  from  four  point  estimates  of  the  function  Y  . 

Each  point  estimate  can  be  viewed  as  a  deterministic  calculation  of  the 
dependent  randan  variable  Y  .  The  equations  presented  above  can  be 
generalized  to  n  random  variables  requiring  2n  point  estimates. 
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CHAPTER  3 


APPLICATION  OF  THE  ANALYSIS  TECHNIQUE 

3.1  MASS-VOLUME  AND  INDEX  DATA 

Dry  densit”  data  obtained  for  an  area  located  near  Yuma,  Arizona 
(Area  A)  are  plotted  versus  depth  in  Figure  3.1-  Four  layers  were  used 
as  an  initial  attempt  to  determine  a  subsurface  profile  based  on  the  dry 
density  data.  The  depth  to  layer  interfaces  and  dry  density  mean  and 
standard  deviation  values  for  each  layer  are  in  Table  3.1.  The  calcu¬ 
lated  F  values  (Equations  2. 1-2. 3)  and  values  for  an  a  of  0.05  for 
the  F  distribution  (  Fq  q,-  )  are  also  shown.  The  one-way  analysis  of 
variance  tests  between  adjacent  layers  is  conducted  by  comparing  the 
calculated  F  value  with  the  value.  The  results  for  the  tests 

U.  UD 

between  layers  1  and  2  and  between  layers  3  and  H  (Table  3.1)  indicate 
that  the  difference  between  the  mean  values  is  not  statistically  signif¬ 
icant,  i.e.,  F1-2  and  F^  are  less  than  their  respective  values  of 

Frt  .  An  attempt  that  uses  a  slightly  different  interface  between 
U .  Ub 

layers  1  through  3  and  combines  layers  3  and  H  is  summarized  in 
Table  3-2. 

In  Table  3.2,  the  test  between  layers  1  and  2  does  not  indicate  a 
statistically  significant  difference,  but  the  test  between  layers  2  and 
3  does  indicate  a  statistically  significant  difference.  Therefore, 
layers  1  and  2  were  combined  and  the  test  repeated.  The  results  of  this 
analysis  are  shown  in  Table  3.3.  Since  the  analysis  of  variance  test 
between  the  two  layers  indicates  a  statistically  significant  difference, 
this  sequence  of  layers  is  used  to  describe  the  dry  density-based 
profile  at  Area  A.  The  results  of  the  analysis  are  shown  graphically  in 
Figure  3.1.  The  results  from  applying  the  procedure  to  the  dry  density 
data  obtained  from  an  area  that  is  adjacent  to  Area  A  (designated 
Area  B)  are  summarized  in  Table  3.^  and  shown  graphically  in  Figure  3.2. 

No  significant  trends  could  be  identified  for  water  content  and 
specific  gravity  versus  depth.  Mean  and  standard  deviation  values  for 
water  content  at  Area  A  are  3.8  and  1.6  percent,  respectively;  for 
Area  B,  these  values  are  3* 5  and  1.3  percent,  respectively.  Mean  and 
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standard  deviation  values  for  specific  gravity  of  2.68  and  0.01,  respec¬ 
tively,  are  used  for  both  areas. 

The  mean  and  standard  deviation  values  for  wet  density  Y  ,  degree 
of  saturation  S  ,  and  air  voids  content  Va  are  calculated  using  the 
Rosenbiueth  procedure.  A  symmetrical  distribution  and  no  correlation  is 
assumed.  From  a  visual  inspection  of  Figures  3.1  and  3.2,  the  assump¬ 
tion  of  a  symmetrical  distribution  for  dry  density  is  probably  adequate 
since  an  equal  amount  of  data  is  approximately  on  either  side  of  the 
mean  value.  Since  the  layering  that  was  determined  for  the  dry  density 
data  could  not  be  applied  to  the  water  content  and  specific  gravity 
data,  it  is  reasonable  to  assume  that  these  properties  are  uncorrelated; 
therefore,  the  value  of  p  in  equations  2.13  and  2.15  is  assumed  to  be 
zero.  The  value  of  Y  in  terms  of  dry  density  Y^  and  water  content 
a)  ,  and  the  value  of  S  and  Va  in  terms  of  Yrf  ,  to  ,  and  specific 
gravity  G„  are  given  by  (Reference  18) 


s  = 


'3.2) 


V 

a 


=  1 


(3.3) 


Equations  3.1-3- 3  can  be  easily  incorporated  into  the  Rosenbiueth 
procedure.  The  resulting  set  of  equations  for  E[Y]  and  V[ Y]  are 

E[Y]  =  J-  [  Y++  +  Y+~  +  Y_v  +  Y“]  ( 3 .*0 


1  ++  2 
V[Y]  =  \  [  ( Y  ) 


( Y+_ ) 2  >  (Y~V  +  (Y_V  ]  -  E2[Y] 


(3.5) 


The  resulting  set  of  equations  for  E[S]  ,  V[S]  ,  EC V_ ]  ,  and  V[V„]  are 

cl  ci 
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ECS]  =  i  [s++++  s++_  *  s+'+  +  s~++  +  s'"  +  s"'+ 


—  +  —  + - .. 

+  s  +  S  J 


(3.6) 


1  +++  2  ++“  2  +— +  2  — ++  2 

V[S]  =  g-  CCS  )  +  (S  )  +  (S  )  +  (S  )  - 

(S  A  (s~V  +  (s"+“)2  +  (S+")2]  -  E2  CS] 


(3.7) 


1  _  +++  ++—  + — +  — ++  — - 

ECv  ]  =  4-  Cv  +  v  +  v  +  v  +  v 

a  8  a  a  a  a  a 


v'+~  +  V+  ] 
a  a 


--+ 

+  V  + 
a 


(3.8) 


1  r, ..++  +  , 


+-+ , *  ,  -  +  +  , 


VCV  ■  8  C(Va  >  +  (Va  >  *  (Va  >  *  (Va  >  * 


(V  )%  (V'  V*  (v/V*  (V*-)2]  -  E2  [VJ 
a  a  a  a  a 


(3.9) 


In  equations  3- ^“3. 9,  terms  with  a  K  +  "  in  the  superscript  indicate  the 

A 

use  of  the  mean  +  1  standard  deviation  (a)  value  and  terms  with  a  "  -  " 
in  the  superscript  indicate  the  use  of  the  mean  -  a  value.  The  proce¬ 
dure  is  illustrated  in  Figure  3-3  for  the  first  layer  of  Area  A.  'f'he 
mean  and  standard  deviation  values  for  the  mass-volume  relations  of  both 
areas  are  summarized  in  Tables  3-5  and  3.6. 

3.2  MECHANICAL  PROPERTY  DATA 

The  mechanical  property  data  associated  with  the  first  layer  of  the 
dry  density-based  profile  for  Area  A  will  be  U3ed  to  illustrate  the 
technique  for  developing  probabilistic  mechanical  properties. 

3.2.1  Grouping  of  the  Mechanical  Property  Data 

The  results  from  uniaxial  strain  (UX),  UX  with  lateral  stress 
measurement  (Ko),  and  hydrostatic  compression  (HC)  tests  conducted  on 
specimens  obtained  f^om  the  upper  49  ft  of  Area  A  were  grouped  to  gener¬ 
ate  the  UX,  Ko,  and  HC  data  bases.  The  results  from  these  tests 
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conducted  on  specimens  obtained  from  the  upper  49  ft  of  Area  B  and  whose 
dry  density  is  within  one  standard  deviation  of  the  mean  dry  density 
determined  in  Section  3.1  were  used  to  supplement  the  data  bases. 
Insufficient  data  are  available  to  group  the  results  from  the  triaxial 
compression  (TXC)  tests  with  respect  to  depth.  The  available  TXC  data 
were  grouped  based  on  the  confining  pressure  at  the  start  of  the  shear 
phase.  Pertinent  information  on  the  individual  tests  is  summarized  in 
Tables  3-7-3.12. 

The  Ko  tests  are  separated  into  two  groups,  i.e.,  Ko  and  EXT  Ko 
tests  (Tables  3-8  and  3.9,  respectively).  The  Ko  tests  were  conducted 
in  a  standard  TX  test  device,  and  the  EXT  Ko  tests  were  conducted  in  an 
extension  TX  test  device  (Reference  19).  The  extension  TX  test  device 
differs  from  the  standard  TX  test  device  in  that  the  piston  is  the  same 
diameter  as  the  test  specimen.  During  an  EXT  Ko  test,  the  axial  stress 
can  be  unloaded  to  zero. 

The  results  from  static  and  dynamic  tests  are  used  in  the  applica¬ 
tion  of  the  analysis  technique.  Static  tests  are  conducted  with  times 
to  peak  stress  of  several  minutes.  Dynamic  tests  are  conducted  with 
times  to  peak  stress  of  several  milliseconds. 

3.2.2  Smoothing  of  the  Data 

The  mechanical  property  test  data  were  smoothed  using  the  technique 
described  in  Section  2.2.1.  The  smoothing  process  begins  by  determining 
which  pare  of  the  data  requires  smoothing.  It  is  suggested  that  a  small 
value  for  n  be  used  initially.  The  section  can  be  smoothed  repeatedly 
with  the  same  window  size  until  a  satisfactory  response  is  produced  or 
until  further  repetitions  produce  no  change  in  the  "smoothed"  curve.  If 
additional  smoothing  is  required,  a  larger  window  size  must  be  used. 

Much  of  the  UX  and  higher  pressure  Ko  and  TXC  test  results  required 
some  smoothing.  Figure  3.4  shows  a  comparison  between  raw  data  and 
smoothed  data.  Note  that  it  is  difficult  to  distinguish  individual 
responses  in  the  raw  data  (Figure  3.4a).  Figure  3.4b  shows  the  data 
after  smoothing.  It  is  much  easier  to  see  the  trends  in  the  data  and  to 
distinguish  individual  responses.  The  general  characteristics  of  the 
data  are  not  altered. 
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3-2.3  Covariance  Analysis 

A  program  used  to  conduct  the  covariance  (COV)  analysis  v/as  written 
for  a  Hewlett-Packard  9816  with  approximately  300  kbytes  of  available 
memory  after  the  ^>0  Basic  operating  system  was  loaded.  This  restricted 
the  number  of  increments  that  could  be  used  in  the  analysis.  Also,  the 
run  time  for  the  program  becomes  unrealistically  high  if  a  large  number 
of  increments  is  used.  It  is  desirable  to  use  as  many  increments  as 
possible  since  the  accuracy  of  the  analysis  increases  with  decreasing 
increment  size.  During  the  development  of  the  program,  it  was 
determined  that  approximately  70  increments  provided  good  results  in  a 
reasonable  length  of  time,  i.e.,  a  few  minutes  of  calculation  time  per 
data  set  depending  on  the  number  of  tests. 

The  curves  used  in  the  COV  analysis  must  be  monotonic  in  the  incre¬ 
mented  axis  direction.  The  analysis  of  loading,  unloading,  and  reload¬ 
ing  responses  are  conducted  separately.  When  the  analysis  is  complete, 
the  probabilistic  loading,  unloading,  and  reloading  responses  are  pieced 
back  together. 

The  COV  analysis  can  be  applied  to  test  data  plotted  in  a  form  that 
is  conducive  to  constitutive  modeling,  i.e.,  axial  strain,  volumetric 
strain,  principal  strain  difference,  axial  stress,  principal  stress 
difference,  and  mean  normal  stress,  or  as  plots  of  "measured  data," 
i.e.,  stresses,  loads,  and  deflections.  If  the  COV  analysis  is  con¬ 
ducted  on  the  "measured  data,"  the  mean  and  standard  deviation  responses 
are  converted  into  the  conducive  forms  using  the  Rosenblueth 
procedure.  Both  methods  will  be  applied  and  their  results  compared. 

3.2.3. 1  UX  test  results.  The  results  of  UX  tests  are  usually 
plotted  as  axial  stress  versus  axial  strain.  The  slope  of  the  curve, 
i.e,,  axial  stress/axial  strain,  is  the  constrained  modulus  M  .  Axial 
stress  is  directly  measured  during  the  test  and  axial  strain  ez  is 
calculated  as 


e 
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where  AH  is  the  measured  deflection  and  HQ  is  the  original  height  of 
the  test  specimen.  This  plot  is  in  the  form  that  is  conducive  to 
constitutive  modeling  and  also  reflects  the  "measured  data." 

Only  the  dynamic  UX  test  aata  will  be  analyzed  since  there  is  an 
insufficient  number  of  static  test  data  to  conduct  a  meaningful  statis¬ 
tical  analysis.  A  composite  plot  of  the  loading  dynamic  UX  test  data 
used  in  the  analysis  is  shown  in  Figure  3.5.  Note  that  all  of  the  data 
does  not  extend  to  the  same  stress  level.  To  utilize  as  much  of  the 
data  as  possible,  the  COV  analysis  was  conducted  on  the  four  groups 
indicated  by  the  dashed  horizontal  lines.  All  of  the  data  is  used  in 
the  first  group.  For  the  remaining  groups,  only  those  tests  which  have 
peak  stresses  at  or  above  the  stress  level  indicate!  by  the  horizontal 
lines  are  used.  An  additional  advantage  to  conducting  the  analysis  of 
the  UX  data  i  several  groups  is  that  smaller  increments  are  used  in  the 
lower  stress  region.  This  provides  a  more  accurate  analysis  in  the 
region  in  which  the  most  curvature  is  occurring. 

In  Section  2.2.2,  it  was  noted  that  the  COV  analysis  could  be 
conducted  by  incrementing  along  the  X-axis  or  the  Y-axis,  i.e.,  the 
strain-  or  stress-axis.  For  the  UX  test,  the  stress-axis  is  chosen 
since  a  stress  is  applied  to  the  specimen,  and  tne  resulting  deflection 
is  measured.  The  applied  stress  is  an  independent  random  variable,  and 
the  resulting  deflection  is  a  dependent  random  variable.  Instead  of 
using  the  incremental  M  ,  compliance  1 /M  is  used  in  the  COV  analysis. 

The  groups  of  loading  dynamic  UX  test  data  used  in  the  COV  analysis 

A 

are  shown  in  Figure  3.6.  The  mean  and  mean  ±  o  responses  are  shown  as 
thick  solid  and  dashed  lines,  respectively.  If  the  data  are  assumed  to 
be  normally  distributed,  then  the  probability  that  a  test  response  is 

A  A 

within  3o  of  the  mean  is  approximately  100  percent.  The  mean  +  3o  re- 

/V 

sponse  may  be  valid;  however,  the  mean  -  3o  response  for  low  stresses  is 
invalid  because  it  leads  to  negative  axial  strains,  which  are  not  possi¬ 
ble  during  a  UX  test.  This  implies  that  the  low  stress  dynamic  UX  data 
are  not  normally  distributed. 

The  mean  and  mean  ±  o  responses  from  Figures  3.6a~3.6b  are  also 
plotted  in  Figure  3.'7a.  Although  the  number  of  curves  decreases  with 
increasing  stress,  the  general  shapes  of  the  probabilistic  responses  do 
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not  change.  A  continuous  probabilistic  response  constructed  with  the 
individual  responses  in  Figure  3.7a  is  shown  in  Figure  3.7b.  The  com¬ 
posite  response  is  constructed  by  using  the  probabilistic  response  for 
group  1  as  the  base.  The  next  section  is  drawn  by  shifting  the  probabi¬ 
listic  response  for  group  2  oetween  axial  stresses  of  12  MPa  and  30  MPa 
to  line  up  with  the  lower  pressure  response.  The  process  is  repeated 
for  each  successive  section.  This  method  of  constructing  the  probabi¬ 
listic  response  is  valid  for  three  reasons:  (1)  by  shifting  the  re¬ 
sponses  in  this  manner,  the  slope  of  the  individual  pieces  is  not  al¬ 
tered;  (2)  as  stress  increases,  the  appearance  of  the  individual  re¬ 
sponses  becomes  very  similar  (Figure  3.5);  and  (3)  the  effect  of  the 
high  stress  tests  on  the  low  stress  probabilistic  responses  is  included 
in  the  analysis,  A  more  rigorous  statistical  analysis  is,  thereby, 
conducted  in  the  region  where  the  greatest  change  in  slope  is  occurring. 

A  composite  of  the  two  sets  of  available  unloading  data  for  the 

A 

dynamic  UX  tests  along  with  the  mean  and  mean  ±  a  responses  are  shown  in 
Figure  3-8.  The  test  data  for  groups  U1  and  U2  are  arbitrarily  normal¬ 
ized  to  axial  strains  of  H  and  8  percent,  respectively.  Unloadings  that 
begin  at  approximately  the  same  stress  level  are  grouped  together  to 
form  the  data  set.  Since  the  unloadings  in  a  given  group  do  not  begin 
at  exactly  the  same  stress  level,  it  is  necessary  to  limit  the  starting 
stress  for  the  unloadings  to  the  lowest  stress  level  for  a  group.  The 
loops  in  the  data  are  caused  by  the  unload-reload  effect  generated  by 
reflected  stress  waves  within  the  test  chamber.  This  caused  the  unusual 
probabilistic  response  near  the  end  of  unloading  U2.  It  is  not  unusual 
to  see  an  unloading  "hook"  in  the  UX  compressibility  response.  However, 
it  is  unusual  to  see  a  response  in  which  the  strain  decreases,  in¬ 
creases,  and  then  decreases  again  near  the  end  of  an  unloading.  For 
this  reason,  the  responses  calculated  for  set  U2  of  the  unloading  data 
vail  be  used  as  only  a  guide.  Probabilistic  reloading  responses  are  not 
provided  due  to  insufficient  data. 

The  constructed  load-unload  probabilistic  dynamic  UX  response  is 
shown  in  Figure  3-9.  A  problem  in  the  construction  process  is  where  to 
place  the  unloading  responses.  From  a  visual  inspection  of  the  test 
data,  the  stiffer  unloadings  tend  to  coincide  witn  tne  surfer  ioaaings, 
and  softer  unloadings  tend  to  coincide  with  softer  loadings. 


The  unloading  responses  are  placed  on  the  appropriate  loading  response 
at  the  starting  stress  level  for  the  unloading  group.  For  example, 
since  the  first  unloading  begins  at  a  vertical  stress  of  29.5  MPa,  the 
mean  unloading  curve  is  shifted  to  begin  at  29.5  MPa  on  the  mean  loading 
curve.  The  stiffest  unloading  curve  is  shifted  to  start  at  29.5  MPa  on 
the  stiffest  loading  curve,  and  the  softest  unloading  curve  begins  at 
29.5  MPa  on  the  softest  loading  curve. 

3. 2. 3. 2  Ko  test  results.  The  Ko  test  data  were  analyzed  by 
(a)  conducting  the  COV  analysis  on  data  that  is  plotted  in  a  form  conve¬ 
nient  for  constitutive  modeling,  i.e.,  axial  stress  versus  axial  strain 
and  principal  stress  difference  (PSD)  versus  mean  normal  stress  (MNS) 
and  (b)  conducting  the  COV  analysis  on  plots  of  data  that  are  closer  to 
that  actually  measured  during  the  Ko  test.  In  Section  3.2.^,  the  re¬ 
sults  from  the  second  analysis  are  transformed  into  the  convenient  forms 
for  constitutive  modeling  using  the  Rosenblueth  procedure.  For  standard 
Ko  tests,  the  "measured  data"  plots  are  load/original  area  (L/Ao)  versus 
axial  strain  and  L/Ao  versus  radial  stress.  The  EXT  Ko  tests  are  con¬ 
ducted  in  a  test  device  that  has  upper  and  lower  cr.ambers  for  applying 
pressure  (denoted  here  as  UC  and  LC,  respectively).  The  "measured  data" 
plots  are  (UC-LC)*Ap/Ao  versus  axial  strain  and  (UC-LC)*Ap/Ao  versus 
LC.  The  area  of  the  piston,  Ap  ,  is  included  to  account  for  an  upward 
stress  on  the  piston  if  the  specimen  has  a  smaller  diameter  than  the 
piston,  and  a  downward  stress  on  the  specimen  if  the  specimen  has  a 
larger  diameter  than  the  piston.  For  Ko  tests,  the  terms  L/Ao  and  (UC- 
LC)*Ap/Ao  are  equivalent  to  PSD  since  radial  strain  is  essentially  zero. 

The  loading  portions  of  Ko  and  EXT  Ko  test  data  plotted  as  de¬ 
scribed  above  are  shown  in  Figures  3-10-3-13.  The  Ko  test  results  are 
analyzed  in  four  groups  (Figures  3.10  and  3-11).  The  test  results 
pointed  out  with  the  letter  "A"  are  only  used  in  group  ,\.  These  results 
were  developed  from  hand- processed  records  and  do  not  have  the  accuracy 
of  the  other  Ko  tests  at  lower  stress  levels.  The  remaining  test  re¬ 
sults  in  group  H  and  the  test  results  in  group  3  are  not  used  in 
groups  1  and  2  because  they  lack  the  accuracy  of  the  results  at  lower 
stress  levels.  One  group  of  loading  EXT  Ko  test  data  are  analyzed 
(Figures  3.12  and  3-13).  Many  of  the  EXT  Ko  test  results  have  "noisy" 
responses  caused  by  (a)  the  electronics  of  the  measurement  system  and 
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(b)  adjustment  of  the  applied  stress  to  maintain  zero  radial  strain. 

For  this  reason,  the  "noise"  in  the  test  results  was  not  completely 
removed  in  the  smooching  process. 

As  discussed  in  Section  3. 2. 3.1,  the  stress-axis  is  incremented 
when  the  data  is  plotted  as  stress  versus  strain.  In  general,  as  MNS 
increases  the  plot  of  PSD  versus  MNS  (or  stress  path)  begins  to  break 
over.  If  this  response  is  assumed  to  be  typical,  more  of  the  loading 
data  is  used  if  the  MNS-axis  is  incremented.  The  L/Ao-  and  (UC- 
LC)*Ap/Ao-axes  are  incremented  to  provide  a  common  axis  for  applying  the 
Rosenbiueth  procedure. 

A 

The  test  data  (thin  solid  line)  and  mean,  mean  ±  c  ,  and 

A 

mean  ±  3o  responses  from  the  COV  analysis  (thick  solid,  thick  dashed, 
and  thin  dashed  lines,  respectively)  are  presented  in  Figures 

A 

3.12-3.17.  All  of  the  test  data  are  within  the  mean  ±  3a  responses. 
However,  like  the  UX  tests,  the  axes  form  a  boundary  for  the  test 
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results.  Since  the  mean  -  3a  responses  cross  into  the  negative  regions 
at  low  stress  levels,  the  data  are  not  normally  distributed  at  low 
stresses.  The  probabilistic  responses  from  the  analysis  of  the  EXT  Ko 
test  data  (Figures  3.12  and  3.13)  are  much  smoother  than  the  data; 
therefore,  the  "noisy"  plots  did  not  significantly  effect  the  results 
from  the  COV  analysis. 

Composite  plots  of  the  results  from  the  analysis  of  the  Ko  test 
data  are  presented  in  Figures  3.18  and  3.19.  The  results  from  the 
analysis  conducted  on  the  stress  versus  strain  plots  for  groups  1  and  Z 
are  much  softer  than  the  results  for  groups  3  and  4.  This  difference  is 
caused  by  (a)  the  data  for  groups  3  and  H  being  stiffer  than  the  data 
for  groups  1  and  2  and  (b)  the  exclusion  of  the  data  in  groups  3  and  4 
in  the  analysis  of  groups  1  and  2.  The  shapes  of  che  respective  stress 
versus  strain  responses  are  similar.  The  results  from  the  analysis 
conducted  on  the  respective  stress  versus  stress  plots  are  similar  for 
all  groups.  The  continuous  responses  in  Figures  3.18  and  3.19  are 
constructed  as  described  in  Section  3.2.3. 1. 

The  results  from  the  analysis  conducted  on  five  sets  of  unloading 
and  two  sets  of  reloading  Ko  test  data  and  two  sets  of  unloading  and  one 
set  of  reloading  EXT  Ko  test  data  are  shown  in  Figures  3.20-3.27.  Each 
set  of  data  is  normalized  to  an  arbitrary  X-axis  value.  The  PSD-axis 
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was  incremented  in  the  analysis  of  the  Ko  test  data  plotted  as  PSD 
versus  MNS  (Figure  3-20)  to  utilize  more  of  the  data.  The  unloading  Ko 
data  sets  U3,  U4,  and  U5  and  reloading  set  R2  (Figures  3*20-3.23)  con¬ 
sisted  of  four  or  less  curves;  therefore,  the  results  of  the  analyses 
will  be  used  as  only  a  guide. 

Constructed  probabilistic  load-unload  and  load-unload-reload 
responses  are  presented  in  Figures  3*28  and  3.29  for  Ko  tests.  Load- 
unload-rcload  responses  from  the  analysis  of  the  "measured  data"  plots 
are  not  presented  since  they  are  generally  not  used  in  constitutive 
modeling.  The  method  outlined  in  Section  3. 2. 3.1  wa3  used  to  place  the 
Ko  unloading  responses  on  the  loading  responses  (Figure  3-28).  For 
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clarity,  the  mean  and  mean  ±  a  responses  for  Ko  unloading  groups  02  ana 
U4  are  plotted  as  thick  solid  and  dashed  lines.  Different  axes  were 
incremented  for  the  loading  and  unloading  test  data  plotted  as  PSD 
versus  MNS.  The  method  used  to  reconstruct  the  load-unload  PSD  versus 

A 

MNS  responses  can  place  the  mean  ±  o  unloadings  an  unequal  distance  from 
the  mean  unloading.  The  unloading  responses  from  group  U4  begin  at 
stress  levels  that  are  higher  than  the  peak  of  the  loading  responses. 
Therefore,  the  starting  locations  for  the  unloadings  are  determined  by 
extrapolating  the  loading  responses  to  the  required  stress  level.  The 
starting  values  for  the  Ko  reloadings  (Figure  3. 29)  are  the  ending 
values  for  the  appropriate  unloadings. 

Each  EXT  Ko  test  was  a  complete  load-unload-relcad  test.  There¬ 
fore,  the  beginning  values  for  the  probabilistic  unloadings  and  reload¬ 
ings  were  the  mean  and  standard  deviation  values  for  the  non- incremented 
axis  at  the  beginning  value  of  the  incremented  axis  (Figure  3.30).  For 
example,  the  mean  and  standard  deviation  values  for  axial  strain  at  the 
beginning  axial  stress  of  unloading  group  U1  (17.3  MPa)  are  9.01  and 
2.79  percent,  respectively.  The  mean  axial  stress  versus  axial  strain 
response  for  group  U1  begins  at  an  axial  strain  of  9.01  percent,  the 
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mean  +  o  response  begins  at  6.22  percent,  and  the  mean  -  o  response 
begins  at  11.80  percent. 

As  stated  earlier,  when  the  Ko  and  EXT  Ko  test  data  are  plotted  as 
PSD  versus  MNS,  more  of  the  data  is  used  if  the  MNS-axis  is 
incremented.  The  effect  of  incrementing  along  the  MNS-axis  as  opposed 
to  the  PSD-axis  is  illustrated  in  Figure  3.31  for  groups  1-4  of  the  Ko 
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test  results.  It  is  clear  from  this  figure  that  incrementing  along  the 
MNS-axis  provides  essentially  the  same  results  as  incrementing  along  the 
PSD-axis. 

3. 2. 3. 3  KC  test  results.  The  results  from  static  and  dynamic 
hydrostatic  compression  (HC)  tests  will  now  be  analyzed.  The  static 
tests  have  times  to  peak  stress  of  several  minutes  and  the  dynamic  tests 
have  times  to  peak  stress  of  approximately  100  milliseconds.  Host  of 
the  data  are  from  the  HC  phases  of  triaxial  compression  tests.  For 
partially  saturated  specimens,  the  volumetric  strain  is  calculated  using 
the  axial  and  radial  strains.  Several  methods  are  available  for  calcu¬ 
lating  the  volumetric  strain  based  on  an  assumed  shape  of  the  specimen 
during  the  test  (Reference  20).  In  general,  the  volumetric  strains 
during  a  HC  test  are  calculated  assuming  deformation  as  a  right  circular 
cylinder  (  dV/Vo(cyi)  ). 

For  HC  test  data,  the  plot  that  is  convenient  for  constitutive 
modeling  is  MNS  versus  dV/Vo(cyl)  where  MNS  is  equivalent  to  the  confin¬ 
ing  pressure.  The  "measured  data"  are  plotted  as  confining  pressure 
versus  axial  strain  and  confining  pressure  versus  radial  strain.  As 
discussed  in  Section  3.2.3. 1,  the  stress-axis  is  incremented. 

The  loading  portions  of  the  static  and  dynamic  test  data  are 
plotted  as  described  above  in  Figures  3-32-3.35.  The  static  test  data 
are  analyzed  in  four  groups  (Figures  3-32  and  3-33).  The  data  in  groups 
three  and  four  are  stiffer  than  the  data  in  the  lower  stress  levels. 
However,  in  the  upper  stress  regions,  there  is  little  change  in 
curvature.  Continuous  responses  constructed  from  the  results  of  the 
analyses  of  each  group  should  not  be  significantly  different  from  the 
results  of  an  analysis  in  which  all  of  the  data  extended  to  the  stress 
level  of  group  4.  The  dynamic  test  data  are  analyzed  in  two  groups. 

The  results  from  the  analysis  conducted  on  the  data  in  the  second  group 
will  have  no  statistical  significance  since  there  are  only  two  curves  in 
the  group. 

The  test  data  (thin  solid  line)  and  mean,  mean  ±  o  ,  and 
mean  ±  3o  responses  from  the  COV  analysis  (thick  solid,  thick  dashed, 
and  thin  dashed  lines,  respectively)  of  the  loading  portions  of  the  HC 
data  are  presented  in  Figures  3.36-3.41.  All  of  the  data  are  within 
38  of  the  mean  response.  Like  the  UX  tests,  the  axes  form  a  boundary 
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for  the  test  results.  Since  the  saean  -  3 <J  responses  e.  oss  into  the 
negative  region  at  low  stress  levels,  the  data  are  not  normally  distrib¬ 
uted  at  low  stresses. 

Composite  plots  of  the  mean  and  mean  ±  o  responses  from  the  analy¬ 
sis  cf  the  loading  static  HC  test  data  are  presented  in  Figures  3-ft2  and 
3.ft3.  The  probabilistic  response  for  group  1  is  softer  than  those  for 
the  other  groups,  and  the  probabilistic  response  for  group  2  is  slightly 
softer  than  that  for  groups  3  and  ft.  The  HC  data  plotted  in 
Figures  3-32  and  3-33  also  indicate  this  order  of  stiffness.  The  con¬ 
tinuous  responses  in  Figures  3*ft2,  3-ft3b,  and  3-ft3d  sre  constructed  by 
the  method  outlined  in  Section  3-2. 3-1-  From  visual  inspection,  the 
responses  shown  in  Figures  3-ftO,  3-ftl.  and  3-ft3  indicate  that  the  mate¬ 
rial  responds  isotropically  to  HC  loading. 

The  data  and  results  from  the  analysis  of  two  groups  of  unloading 
and  reloading  static  HC  test  data  and  two  groups  of  unloading  dynamic  HC 
test  data  are  presented  in  Figures  3-ftft_3-ft8.  The  responses  are  normal¬ 
ized  to  arbitrary  values  of  strain.  There  are  more  unloading  responses 
in  group  U2  of  the  dynamic  HC  tests  than  there  are  in  group  2  of  the 
loading  responses.  Group  U2  consists  of  all  available  dynamic  HC 
unloadings.  The  groups  of  loading  dynamic  HC  test  data  included  only 
the  tests  with  times  to  peak  stress  of  approximately  100  msec.  The 
loading  portion  of  the  faster  dynamic  HC  tests  were  very  erratic,  but 
the  unloadings  were  well  behaved. 

Constructed  load- unload -reload  responses  are  shown  in  Figures  3-ft9 
and  3-50.  The  procedure  ^scribed  in  Section  3- 2.3-2  was  used  to  deter¬ 
mine  where  the  unloadings  and  reloadings  should  be  placed.  Although 
there  are  only  two  curves  in  group  2  of  the  loading  dynamic  HC  test 
data,  the  constructed  response  looks  reasonable  (Figure  3-50). 
Constructed  load- unload-reload  responses  for  the  plots  of  "measured 
data"  are  not  presented  since  they  are  not  generally  used  in  constitu¬ 
tive  modeling. 

3. 2. 3- ft  TXC  test  results.  Tne  results  from  TXC  tests  are  plotted 
as  PSD  versus  axial  strain  and  PSD  versus  principal  strain  difference 
for  use  in  constitutive  models.  For  subsequent  use  in  the  Rosenblueth 
technique,  the  "measured  data"  plots  are  L/Ao  versus  axial  strain  and 
L/Ao  versus  radial  strain.  When  plotted  as  PSD  versus  axial  strain  and 
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PSD  versus  principal  strain  difference,  the  strain-axis  is  incremented 
to  utilize  more  of  the  data.  From  Figure  3.51,  incrementing  the  strain- 
axis  as  opposed  to  the  PSD-axis  has  little  effect  on  the  results  of  the 
GOV  analysis.  The  L/Ao-axis  is  incremented  in  the  analysis  of  the 
"measured  data"  plots. 

The  data  and  results  from  the  COV  analysis  conducted  on  the  loading 
portion  of  the  static  TXC  test  data  plotted  as  described  above  are  shown 
in  Figures  3.52-3.6S  for  confining  pressures  ranging  from  zero  to 
80.0  MPa.  The  test  data  are  normalized  to  zero  strain  prior  to  conduct¬ 
ing  the  COV  analysis.  With  the  exception  of  the  data  for  confining 
pressures  of  17.5,  20.0,  and  55.0  MPa  (Figures  3*55,  3*56,  3-58,  3-64, 
3-65,  and  3-67),  there  is  insufficient  data  to  conduct  meaningful,  sta¬ 
tistical  analyses.  The  radial  strains  are  plotted  as  negative  values  in 

Figures  3-61-3-69  because  radial  deflection  during  the  shear  phase  of  a 
TXC  test  is  generally  expansive.  The  radial  strain  from  some  of  the 
tests  are  positive  during  shear.  This  response  is  likely  caused  by 
wrinkles  in  the  membrane  surrounding  the  specimen,  irregularities  in  the 
specimen  surface,  etc.  Most  of  the  test  data  are  within  the  mean  ±  3a 
responses.  The  data  are  not  normally  distributed  at  low  values  of 
strain  since  the  mean  -  3cr  responses  in  Figures  3-55-3-60  and  3-64-3.69 
imply  negative  values.  Some  of  the  individual  tests  had  "noisy"  re¬ 
sponses  that  do  not  significantly  effect  the  calculation  of  the  mean  and 
mean  ±  0  responses. 

The  data  and  results  from  COV  analysis  of  the  available  unloading 
and  reloading  static  TXC  test  data  are  presented  in  Figures  3.70-3-85 
for  the  various  confining  pressures.  For  a  given  confining  pressure, 
the  responses  are  grouped  based  on  the  strain  at  which  the  unloadings 
begin.  The  stress-axis  is  incremented  to  utilize  more  of  the  unloading 
and  reloading  data.  Much  less  variation  occurs  in  the  unloading  and 
reloading  data  than  in  the  leading  data. 

The  constructed  probabilistic  load- unload-reload  responses  for  PSD 
versus  axial  strain  and  PSD  versus  principal  strain  difference  plots  of 
the  TXC  test  data  are  presented  in  Figures  3-86-3.93.  The  beginning 
strains  for  the  unloading  and  reloading  responses  are  the  average  values 
for  the  respective  groups  of  data.  The  unloadings  and  reloadings  do  not 
extend  to  the  loading  responses  since  the  unloading  (or  reloading)  which 
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starts  at  the  lowest  value  of  PSD  (highest  for  reloadings)  governs  the 
starting  PSD  for  the  group.  The  final  unloading  responses  start  at 
values  of  strain  that  are  larger  than  the  end  of  the  loading  response. 
The  loading  responses  must  stop  at  the  smallest  ending  value  of  strain 
for  a  given  group  of  loading  test  data.  The  load- unload- reload 
responses  for  the  "measured  data"  plots  are  not  constructed  since  they 
are  not  used  in  constitutive  modeling. 

Results  from  the  COV  analysis  in  which  axial  strain  reached  15  per¬ 
cent  or  a  "peak"  strength  was  obtained  do  not  necessarily  provide  mean 
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and  mean  ±  c  values  for  strength  at  failure.  For  example,  the  results 
for  a  confining  pressure  of  0.35  MPa  (Figure  3.53)  imply  mean  and 
mean  ±  o  values  for  "peak"  strength  of  1.4*1  MPa,  1.79  MPa,  and  1.34  MPa, 
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respectively.  The  implied  mean  +  a  and  mean  -  o  values  are  not  an  equal 
distance  from  the  mean.  A  probabilistic  "peak"  strength  envelope  based 
on  the  test  data  will  be  provided  in  Chapter  4. 

3.2.4  Rosenblueth  Procedure 

In  this  section,  the  Rosenblueth  equations  presented  in  Chapter  2 
(equations  2.15  through  2.18)  are  used  to  convert  the  results  from  the 
COV  analysis  conducted  on  plots  of  "measured  data"  (Section  3-2.3)  into 
forms  that  are  conducive  to  constitutive  modeling.  The  mean  and  vari¬ 
ance  of  secant  slopes  are  used  in  the  conversion.  Since  slopes  are 
calculated,  the  effect  of  correlation  coefficients  on  the  results  from 
the  Rosenblueth  procedure  must  be  evaluated.  The  calculation  of  the 
correlation  coefficients  is  discussed  in  Section  3-2. 4.2. 

3.2.4. 1  Equations  for  the  function  y(x).  The  equations  for  con¬ 
verting  the  probabilistic  responses  from  the  COV  analysis  conducted  on 
the  Ko,  EXT  Ko,  HC,  and  TXC  "measured"  test  data  are  presented  in 
Figures  3.94-3.97,  respectively.  For  Ko  tests,  the  target  plots  are  PSD 
versus  MNS  (stress  path)  and  axial  stress  versus  axial  strain 
(Figure  3-94).  The  stress  path  is  developed  by  simply  modifying  the 
L/Ao  versus  radial  stress  results.  The  mean  and  variance  of  the  inverse 
slope  are  calculated  (MNS/PSD).  The  probabilistic  response  is 
calculated  by  multiplying  the  resu] ts  of  the  Rosenbluetn  procedure  by 
L/Ao  since  L/Ao  is  equal  to  PSD  for  a  Ko  test.  The  Rosenblueth 
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procedure  requires  two  combinations  since  there  is  only  one  dependent 
variable  (radial  stress).  The  axial  stress  versus  axial  strain  response 
requires  the  COV  analysis  results  for  data  plotted  as  L/Ao  versus  axial 
strain  and  L/Ao  versus  radial  stress.  There  are  four  combinations  since 
there  are  two  dependent  variables.  The  mean  and  variance  of  the  inverse 
slope  are  calculated  (axial  strain/axial  stress).  The  probabilistic 
values  for  axial  strain  are  calculated  for  values  of  axial  stress. 

Axial  stress  is  calculated  using  the  given  value  for  L/Ao  and  the  mean 
value  for  radial  stress.  The  reason  for  using  the  inverse  slope  is  the 
same  as  that  described  in  Section  3.2. 3. 2  a  or  CO V  analysis  on  data 
plotted  as  axial  stress  versus  axial  strain,  i.e.,  strain  is  the 
dependent  variable,  and  stress  is  the  independent  variable. 

The  target  plots  for  the  EXT  Ko  tests  are  the  same  as  those  for  the 
Ko  tests,  and  the  procedures  are  the  same  (Figure  3.95).  The  only 
difference  is  the  substitution  of  (UC-LC)*Ap/Ao  for  L/Ao  and  LC  for 
radial  stress. 

The  target  plot  for  the  HC  tests  is  MNS  versus  dV/Vo(cyl).  The 
Rosenblueth  equation  is  developed  in  Figure  3-96.  Mean  normai  stress  i3 
equal  to  the  confining  pressure  since  it  is  an  equal  all-around 
stress.  The  inverse  slope  (  dV/Vo(cyl)  /  MNS  )  is  calculated  since 
stress  is  the  independent  variable  and  strain  is  the  dependent 
variable.  The  mean  and  variance  of  dV/Vo(cyl)  is  calculated  by  multi¬ 
plying  the  mean  and  variance  of  the  inverse  slope  by  the  current  value 
of  confining  pressure  (MNS).  The  Rosenblueth  procedure  requires  four 
combinations  since  there  are  two  dependent  variables  (axial  strain  and 
radial  strain). 

The  target  plots  for  the  TXC  tests  are  PSD  versus  axial  strain  and 
PSD  versus  principal  strain  difference.  The  Rosenblueth  equations  are 
developed  in  Figure  3.97.  A  problem  with  the  TXC  tests  is  that  PSD  is 
calculated  as  load  divided  by  the  current  area,  where  the  current  area 
is  the  original  area  multiplied  by  the  quantity  one  minus  the  true 
radial  strain  squared.  The  plots  of  L/Ao  versus  radial  strain  are 
normalized  to  zero;  therefore,  the  true  radial  strain  is  not  incorpo¬ 
rated  in  the  COV  analysis.  The  results  of  the  COV  analysis  on  data 
plotted  as  L/Ao  versus  radial  strain  are  shifted  to  account  for  the 
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normalization.  The  shift  is  determined  by  calculating  the  mean  and 
standard  deviation  for  radial  strain  at  the  start  of  shear  for  the  tests 
at  a  given  confining  pressure.  The  mean  result  from  the  COV  analysis  is 
shifted  by  the  mean  radial  strain  at  the  start  of  shear.  The  mean  +  a 
and  mean  -  c  results  are  shifted  by  the  mean  +  a  and  mean  -  o  values  for 
radial  strain  at  the  start  of  shear,  respectively.  The  inverse  slope  is 
calculated  (strain/PSD)  since  PSD  is  the  independent  variable  and  strain 
is  the  dependent  variable.  The  mean  and  variance  for  strain  are  calcu¬ 
lated  by  multiplying  the  mean  and  variance  of  the  slope  by  the  mean 
value  of  PSD  for  the  current  increment.  The  mean  value  of  PSD  is  calcu¬ 
lated  using  the  current  value  of  L/Ao  and  the  mean  true  radial  strain. 
The  Rosenblueth  procedure  requires  four  combinations  since  there  are  two 
dependent  random  variables  (axial  strain  and  radial  strain). 

3. 2. 4.2  Correlation  coefficients.  In  this  report,  correlation 
coefficients  (p)  are  estimated  by  a  linear  regression  analysis  conducted 
on  the  two  random  variables.  The  shortcoming  of  this  approach  is  th3t 
large  absolute  values  of  p  imply  that  the  degree  of  correlation  be¬ 
tween  the  two  random  variables  is  high,  whereas  small  absolute  values 
imply  only  the  weakness  of  a  linear  trend  and  not  necessarily  a  weakness 
in  correlation.  Positive  values  of  p  tend  to  imply  that  large  values 
of  cne  variable  coincide  with  large  values  of  the  other  variable  and 
visa  versa.  Negative  values  of  p  tend  to  imply  that  large  values  of 
one  variable  coincide  with  small  values  of  the  other  variable  and  visa 
vers? . 

From  equations  2.13  and  2.15,  values  of  p  approaching  1  will 
increase  the  significance  of  the  ++  and  —  calculations,  and  values  of 

p  approaching  -1  will  increase  the  significance  of  the  +-  and 

% 

calculations.  Equal  weight  is  given  to  all  parts  of  the  Rosenblueth 
calculation  when  p  i3  equal  to  zero.  Depending  on  the  function  yix) 
and  the  values  of  the  correlation  coefficients,  the  effect  of  including 

p  in  the  calculations  can  be  significant.  The  effect  of  p  on  the 
mean  response  varies  depending  on  the  relative  magnitudes  of  the  ++,  +-, 
-+,  and  —  values.  The  effect  of  p  on  the  variance  is  predetermined 
by  the  sign  of  p  and  the  rank  of  the  various  combinations.  If  the  ++ 
and  —  combinations  provide  the  extremes,  the  variance  will  increase  for 
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positive  values  of  p  and  decrease  for  negative  values  of  p  .  If  the 
+-  and  -+  combinations  provide  the  extremes,  then  the  variance  will 
decrease  for  positive  values  of  p  and  increase  for  negative  values 
of  p  . 

Given  n  pairs  of  values,  the  correlation  coefficient  is  estimated 
by  the  sample  correlation  coefficient  r 
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Correlation  coefficients  for  the  loading  portion  of  each  test  type 
are  plotted  versus  the  incremented  axis  in  Figures  3.98-3.103.  For  Ko 
tests  (Figure  3.98),  the  correlation  coefficients  for  groups  1  and  2  are 
between  0,5  and  1.0  to  a  L/Ao  of  3.1*  MPa;  the  correlation  coefficients 
are  approximately  0.75  above  a  L/Ao  of  3.4  MPa.  For  groups  3  and  4  of 
the  Ko  test  data,  the  values  for  p  are  initially  between  0.7  and 
1.0.  The  values  of  p  decrease  rapidly  to  0  at  an  axial  stress  of 
14  MPa.  Above  14  MPa,  the  values  for  group  3  fluctuate  about  0  while 
the  values  for  group  4  continue  to  decrease  to  a  low  of  -0.4  at 
66  MPa.  The  correlation  coefficients  for  EXT  Ko  tests  (Figure  3.99) 
gradually  decrease  from  0.95  to  0.4  as  (UC-LC)*Ap/ Ao  increases. 

The  correlation  coefficients  for  groups  1  and  2  of  the  static  HC 
test  data  (Figure  3.100a)  decrease  rapidly  as  confining  pressure  in¬ 
creases  to  1  MPa.  Above  a  confining  pressure  of  1.5  MPa,  the  correla¬ 
tion  coefficients  for  group  1  are  approximately  0.35;  for  group  2,  the 
correlation  coefficients  are  approximately  0  above  6  MPa. 
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The  correlation  coefficients  for  groups  3  and  4  (Figure  3.100b)  change 
rapidly  below  a  confining  pressure  of  4  MPa.  For  group  3,  the  values 
then  decrease  slowly  from  -0.17  to  -0.45  at  32  MPa;  for  group  *1,  the 
correlation  coefficients  are  approximately  0.65  above  4  MPa.  The  plots 
of  confining  pressure  versus  p  for  groups  1  and  2  of  the  dynamic  HC 
tests  are  in  Figure  3.101.  The  correlation  coefficients  for  group  2  are 
a  constant  -1  since  there  are  only  two  tests.  For  group  1 ,  the  value  of 
p  is  approximately  0  initially;  it  then  decreases  to  -0.42  as  confin¬ 
ing  pressure  increases  to  4  MPa.  Above  4  MPa  the  values  of  p  for 
group  1  increase  to  0.55  at  22  MPa. 

Most  of  the  correlation  coefficients  for  the  TXC  tests  decrease  as 
L/Ao  increases  (Figures  3-102  and  3.103).  At  low  values  of  L/Ao,  the 
TXC  tests  at  confining  pressures  less  than  10  MPa  (Figure  3*102)  have 
correlation  coefficients  that  alternate  from  large  positive  to  large 
negative  values.  The  values  for  p  ultimately  approach  -1«  Since 
there  are  only  two  tests  for  a  confining  pressure  of  0.7  MPa,  the  values 
of  p  are  either  1  or  -1.  The  correlation  coefficients  for  TXC  tests 
at  confining  pressures  greater  than  10  MPa  (Figure  3.103)  are  not  as 
erratic  at  low  stresses  as  the  correlation  coefficients  in 
Figure  3.102.  Except  for  the  values  of  p  for  a  confining  pressure  of 
80  MPa  (Figure  3.103d),  the  correlation  coefficients  in  Figure  3.103 
gradually  decrease  as  L/Ao  increases;  the  values  approach  -1  near  peak 
L/Ao. 

3. 2.4. 3  Results.  The  probabilistic  plots  obtained  by  applying  the 
Rosenblueth  procedure  to  the  results  from  the  C0V  analysis  of  the  "mea¬ 
sured  data"  are  shown  in  Figures  3.104-3.120.  For  the  Ko,  EXT  Ko,  and 
HC  tests  (Figures  3.104-3.108),  the  thick  lines  are  the  results  that 
include  correlation  coefficients.  Correlation  coefficients  are  not 
needed  to  calculate  the  stress  path  (Figures  3.105  and  3.106b)  since 
there  is  only  one  variable.  The  inclusion  of  p  had  a  greater  effect 
on  the  variance  than  the  mean  response. 

The  load-unload-reload  responses  calculated  by  applying  the  Rosen¬ 
blueth  procedure,  including  p  ,  to  the  "measured  data"  from  TXC  tests 
are  shown  in  Figures  3.109-3.117  for  the  various  confining  pressures. 

The  unloadings  and  reloadings  are  shifted  to  the  corresponding  mean 
values  of  axial  strain  and  principal  strain  difference  determined  as 
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described  in  Section  3. 2. 3. 4.  As  illustrated  in  Figures  3.118-3.120, 
the  inclusion  of  p  only  affected  the  variance. 

The  probabilistic  responses  from  applying  the  Rosenblueth  procedure 
to  the  constructed  continuous  plots  of  "measured  data"  from  Ko  and  HC 
tests  are  shown  in  Figures  3.121-3.123.  The  value  of  p  for  a  given 
increment  is  determined  by  which  group  that  section  comes  from.  For 
example,  if  a  section  of  the  combined  response  comes  from  group  3,  then 
the  correlation  coefficients  for  that  section  of  the  calculated  response 
also  come  from  group  3.  Where  relatively  severe  shifts  in  p  occur  at 
section  interfaces,  the  inclusion  of  p  results  in  very  "jagged"  re¬ 
sponses;  whereas,  the  omission  of  p  results  in  smooth  responses.  The 
omission  of  p  seems  to  provide  an  "average"  response. 

The  probabilistic  load-unload  and  load- unload-reload  responses  for 
Ko,  EXT  Ko,  and  HC  tests  are  shown  in  Figures  3.12^-3.1 28.  For  the  Ko 
and  HC  tests,  correlation  coefficients  are  included  in  the  calculation 
of  the  unloading  and  reloading  responses,  but  are  not  included  in  the 
calculation  of  the  loading  responses.  The  starting  values  of  axial 
stress  for  the  Ko  unloading  and  reloading  responses  are  calculated  as 
L/Ao  at  the  start  of  the  unloading  or-  reloading  plus  the  mean  radial 
stress  for  all  of  the  unloading  or  reloading  curves  in  a  given  group. 

The  starting  value  of  axial  strain  for  an  unloading  response  is  the 
value  on  the  leading  response  at  the  calculated  value  of  axial  stress. 
The  starting  value  of  MNS  for  the  unloading  stress  path  is  the  inter¬ 
section  of  the  loading  response  and  the  starting  L/Ao  (PSD)  value  fer 
the  group  of  unloading  data.  For  both  plot  types,  a  given  reloading 
response  begins  at  the  end  of  the  associated  unloading  response. 

The  starting  value  of  axial  stress  for  an  EXT  Ko  unload  or  reload 
response  is  calculated  as  the  value  of  (UC-LC)*Ap/Ao  plus  the  mean  LC 
stress  at  the  beginning  of  a  given  group  of  responses.  The  starting 
values  of  axial  strain  are  the  mean  and  mean  ±  o  values  calculated  at 
the  starting  value  of  (UC-LC) *Ap/Ao  for  the  unnormalized  unloading  or 
reloading  group.  The  mean  and  mean  ±  o  values  for  MNS  at  the  beginning 
of  unload  or  reload  PSD  versus  MNS  responses  are ‘cai culated  as  the 
starting  value  of  (UC-LC)* Ap/Ao  divided  by  3  plus  the  mean  and 
mean  ±  a  LC  stress  values  for  the  appropriate  set  of  responses. 
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The  starting  points  for  the  HC  unloading  responses  are  the  MNS  at 
the  start  of  the  unloading  set  and  the  strain  at  which  that  value  of  MNS 
intersects  the  the  loading  responses;  the  reloadings  begin  at  the  end  of 
their  associated  unloading  response.  For  all  test  types,  the  stiff 
unloading  and  reloading  responses  are  placed  with  the  stiff  loading 
response,  and  the  soft  unloading  and  reloading  responses  are  placed  with 
the  soft  loading  response. 

3.2.5  Comparison  of  the  Results  from  the  Two  Analysis  Techniques 

The  results  from  the  COV  analysis  of  the  loading  mechanical  prop¬ 
erty  test  data  plotted  in  forms  conducive  to  constitutive  modeling  (COV 
technique)  and  the  results  from  applying  the  Rosenblueth  procedure  to 
the  results  from  the  COV  analysis  of  ’’measured  data"  (COV/Rosenblueth 
technique)  are  compared  in  Figures  3*1 29-3*  1  41 .  The  thick  lines  are  the 
results  from  the  COV  technique,  and  the  thin  lines  are  the  results  from 
the  COV/Rosenblueth  technique.  Correlation  coefficients  are  not 
included  in  the  Rosenblueth  calculations. 

Both  techniques  provide  similar  mean  responses  when  applied  to  Ko, 
EXT  Ko,  static  HC,  and  dynamic  HC  test  data  (Figures  3-129-3.152, 
respectively).  The  variance  from  the  COV  technique  is  (a)  slightly  less 
than  the  variance  from  the  COV/Rosenblueth  technique  for  Ko  test  data, 
(b)  very  similar  to  the  variance  from  the  COV/Rosenblueth  technique  for 
EXT  Ko  test  data,  (c)  slightly  larger  than  the  variance  from  the 
COV/Rosenblueth  technique  for  static  HC  test  data,  and  (d)  for  dynamic 
HC  test  data  below  a  MNS  of  about  18  MPa,  smaller  than  the  variance  from 
the  COV/Rosenblueth  technique  and,  above  18  MPa,  larger  than  the  vari¬ 
ance  from  the  COV/Rosenblueth  technique. 

For  TXC  tests  at  confining  pressures  less  than  10  MPa  (Figures 
3.133-3.135),  the  responses  from  the  COV/Rosenblueth  technique  have 
softer  mean  responses  and  greater  variances  than  those  from  the  COV 
technique.  For  confining  pressures  greater  than  10  MPa  (Figures 

A 

3. 136-3. 1^1) ,  the  mean  and  mean  ±  0  responses  from  both  techniques  are 
very  similar.  The  greatest  difference  occurs  near  the  PSD  level  where 
the  curves  begin  to  break  over.  From  Section  3-2.3.^,  the  loading  TXC 
test  data  exhibit  one  or  more  of  the  following  phenomena:  (a)  the 
strengths  continue  to  increase  during  the  test  such  tnat  the  shear 


responses  are  concave  to  the  MNS  axis,  (b)  the  strength  for  a  given  test 
begins  to  decrease  after  achieving  a  peak  PSD.  or  (c)  at  some  level  of 
PSD,  the  strength  remains  nearly  constant  with  increasing  strain.  Since 
the  stress  axis  is  incremented  in  the  COV/Rosenblueth  technique,  these 
phenomena  are  not  incorporated  into  the  probabilistic  responses. 

Tne  unloading  and  reloading  responses  from  the  two  techniques  are 
very  similar. 


Table  3*1-  First  attempt  at  the  dry  density-based  profile  for 
Area  A. 


Layer 

No. 

Depth 

ft 

Mean  Dry 
Density 
pcf 

Dry  Density 
Standard 
Deviation 
pcf 

Number  of 
Dry  Density 
Points 

F 

F0.05 

0 

1 

114.85 

5.17 

260 

21 

0.07 

-3.84 

2 

115.00 

4.78 

123 

39 

12.78 

-3.92 

3 

111.09 

5.00 

23 

59 

2.24 

4.13 

4 

108.26 

5.88 

12 
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3.2. 

Second  attempt  at  the  dry  density-based 
Area  A. 

profile  for 

Dry  Density 

Mean  Dry 

Standard 

Number  of 

Layer 

Depth 

Density 

Deviation 

Dry  Density 

No. 

ft 

pcf 

pcf 

Points 

F 

F0.05 

0 

1 

11 4.72 

5.30 

207 

17 

0.07 

-3.84 

2 

114.71 

4.90 

193 

49 

17.34 

-3.92 

3 

109.60 

5.80 

18 

34 


rrrr' 


r  ’rv*Y: 


<-«►*-****•«-  1  r  *  -*  **  ■%; 


»  -y  -n  »  **  w«*v  «  a*  a  ft  a  <  ^  i"-** 


Table  3-3.  Third  attempt  at  the  dry  density-based  profile  for 
Area  A. 


Layer 

No. 

Depth 

ft 

Mean  Dry 
Density 
pcf 

Dry  Density 
Standard 
Deviation 
pcf 

Number  of 
Dry  Density 
Points 

F 

F0.05 

0 

1 

114. 72 

5.11 

400 

A9 

17.09 

-3-92 

2 

109.60 

5.80 

18 

Table  3.  A. 

Dry  density- 

based  profile  for  Area  B. 

Dry  Density 

Mean  Dry 

Standard 

Number  of 

Layer 

Depth 

Density 

Deviation 

Dry  Density 

P 

No. 

ft 

pcf 

pcf 

Points 

F  *0.05 

0 

1 

12 

117.72 

11.79 

78 

30.39  -3.  92 

2 

32 

113.82 

5-  45 

192 

11.38  -3-92 

3 

67 

115.86 

A.  72 

119 

28.59  -3.92 
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Table  3. 10.  Summary  of  the  results  from  static  hydrostatic  compression 
tests  used  in  the  analysis  of  mechanical  properties  for 
layer  1  of  Area  A. 


Depth 

ft  m 

Posttest 

Water 

Content 

a 

p 

Wet 

Density 

g/cc 

Dry 

Density 

g/cc 

Degree  of 
Saturation 
% 

Air 

Voids 

Content 

2 

Confining 

Pressure 

MPa 

4.12 

1.26 

7.62 

1.977 

Area  A 

1.837 

44.5 

17.5 

17.40 

5.40 

1 .65 

3.77 

1 .840 

1.773 

19.8 

27.2 

58.00 

5.50 

1.68 

4.24 

1.770 

1.698 

19.6 

29.4 

20.60 

6.74 

2.05 

8.11 

1.868 

1.728 

39.4 

21 .5 

17.40 

8.70 

2.65 

4.31 

1.934 

1.854 

25.9 

22.8 

17.40 

8.90 

2.71 

4.70 

1.936 

1.849 

28.0 

22.3 

17.40 

9.00 

2.74 

4.09 

1.891 

1.817 

23.1 

24.8 

80.00 

9.20 

2.80 

5.79 

2.041 

1 .929 

39.9 

16.8 

16.90 

10.30 

3.14 

2.83 

1.890 

1.838 

16.6 

26.2 

83.00 

10.50 

3.20 

3.62 

1.994 

1.924 

24.7 

21 .2 

55.00 

11.30 

3-45 

3.82 

1.971 

1.898 

24.9 

21.9 

17.70 

11.40 

3.48 

4.57 

1.937 

1 .852 

27.4 

22.4 

16.40 

13-80 

4.21 

4.69 

1.965 

1.877 

29.4 

21.2 

62.10 

14.03 

4.28 

7.01 

1 .716 
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Table  3-12.  Summary  of  the  results  from  static  trlaxial  compression  tests  used  In  the  analysis 
of  mechanical  properties  for  layer  1  of  Area  A. 


Mass-Volume 

Data 

"Peak" 

Strength 

Posttest 

Air 

Principal 

Mean 

NY  ter 

Wet 

Dry 

Degree  of 

Voids 

Confining 

Stress 

Normal 

Axial 

Depth 

Content 

Density 

Density 

Saturation 

Content 

Pressure 

Difference 

Stress 

Strain 

ft 

ra 

1 

g/cc 

g/cc 

J 

% 

MPa 

MPa 

MPa 

1 

Area  A 
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3.32 

3.65 

1.973 

1.909 

29.0 

22.0 

0.7 

2.9 

1.5 

5.3 

10.10 

3.08 

2.52 

1.873 

1.827 

19.5 

27.2 

0.7 

2.2 

1 .9 

10.9 

1)5. 40 

13.81 

9.66 

1.997 

1.860 

28.3 

21.9 

1.7 

5.  ' 
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95.85 
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1.927 

1.879 

17.6 
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1.7 

5.1 

3-9 

15.0 

55.05 

16.78 

9.77 

1.922 

1.751 

99.3 

17.6 

1.7 

3.3 

2.8 

15.0 

130.60 

39.82 

♦  66 

1.807 

1.727 

22.6 

27.5 

1.7 

9.9 

3.2 

15.0 

11.10 

1.30 

3.11 
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26.8 

17.9 

3.9 

8.8 

6.3 
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3.80 
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1.875 

23.7 
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9.3 

6.6 

12.3 
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27.5 
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6.9 
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9.57 
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1.951 
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C(Yd)  =  114.72  per  ,  o(Yd)  =  5.11  per 
E(w)  =  3-81  %  ,  o(w)  =  1.5 6% 

E(GS)  =  2.68  .  o(Gs)  =  0.01 
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Figure  3.3.  Rosenblueth  procedure  for  calculating  probabilistic  values  \ 

for  y  ,  S  ,  and  from  probabilistic  values  of  ,  ) 

w  ,  and  G 

'  c*  r 


a.  Composite  of  COV  analysis  results. 


Figure  3.7.  Composite  and  constructed  continuous  plots  of  the  results 
from  COV  analysis  of  loading  dynamic  UX  test  data. 
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Figure  3.8.  Results  from  Che  COV  analysis  of  unloading  dynamic  UX  test  data. 
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b.  L/Ao  versus  radial  stress. 

Figure  3.11.  "Measured  data"  plots  :>f  the  loading  portion  of  K  test 
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a.  Composite  plots.  b.  Constructed  continuous  plots. 

Figure  3.19.  Composite  and  constructed  continuous  results  from  the  COV  analysis  of  the  loading  portion 
of  "measured"  K  test  data. 
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a.  Unloading  groups 
Ul,  U2,  and  U3. 


b.  Unloading  groups 
U4  and  U5. 


-  TEST  DRTR 

-  MERN 

-  MERN  ±  (f 


c.  Reloading  groups 
Rl  and  R2. 


Figure  3.20.  Results  from  the  COV  analysis  of  the  unloading  and 

reloading  portions  of  axial  stress  versus  axial  strain 
K  test  data. 
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LOflD/Ro,  MPa 


RXIfiL  STRRIN.  PERCENT 

a.  Unloading  groups 
Ul,  U2 j  and  U3. 


RXIfiL  STRRIN,  PCRCENT 


b.  Unloading  groups 
U4  and  U5. 


-  TEST  DATA 

-  MEAN 

-  MEAN  ±  f 


c.  Reloading  groups 
R1  and  R2. 


Figure  3.22.  Results  from  the  COV  analysis  of  the  unloading  and 

reloading  portions  of  L/Ao  versus  axial  strain  K  test 
da  ta . 
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20.0 


a.  Unloading  groups 
Ul,  U2>  and  U3. 


b.  Unloading  groups 
U4  and  U5. 


c.  Reloading  groups  R1  and  R2. 


-  TEST  DOTH 

-  MEHN 

-  MERN  ±  cr 


Figure  3.23.  Results  from  the  COV  analysis  of  the  unloading  and 

reloading  portions  of  L/Ao  versus  radial  stress  K  test 
data . 
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b.  Reloading  group  Rl. 

Figure  3.24.  Results  from  the  COV  analysis  of  the  unloading  and 

reloading  portions  of  axial  stress  versus  axial  strain 

EXT  K  test  data, 
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20 


a.  Unloading  groups  U1  and  U2. 
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Figure  3-27.  Results  from  the  COV  analysis  of  the  unloading  and 

reloading  portions  of  (UC-LC)*Ap/Ao  versus  LC  EXT  K  test 
data.  ° 
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Figure  3.35.  "Measured  data"  plots  of  the  loading  portion 
IiC  test  data. 
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Figure  3.39.  Results  from  the  COV  analysis  of  the  loading  portion  of 
dynamic  HC  test  data  plotted  in  a  form  that  is  conducive 
to  constitutive  modeling. 
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Figure  3.40.  Results  from  the  COV  analysis  of  the  loading  portion  of 
confining  pressure  versus  axial  strain  dynamic  HC  test 
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Figure  3.41.  Results  from  the  COV  analysis  of  the  loading  portion  of 
confining  pressure  versus  radial  strain  dynamic  liC  test 
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Figure  3.42.  Composite  and  constructed  continuous  results  from  the  COV 
analysis  of  the  loading  portion  of  static  HC  test  data 
plotted  in  a  form  that  is  conducive  to  constitutive 
modeling. 
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Figure  3. A3.  Composite  and  constructed  continuous  results  from  the  COV 
analysis  of  "measured"  static  HC  test  data. 
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Figure  3.44.  Results  from  the  C0V  analysis  of  the  unloading  and 

reloading  portions  of  static  11C  test  data  plotted  in  a 
form  that  is  conducive  to  constitutive  modeling. 
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Figure  3.45.  Results  from  the  COV  analysis  of  the  unloading  and 

reloading  portions  of  confining  pressure  versus  axial 
strain  static  HC  test  data. 
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.46.  Results  from  the  COV  analysis  of  the  unloading  and 
reloading  portions  of  confining  pressure  versus 
radial  strain  static  HC  test  data. 
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Figure  3.47.  Results  from  the  COV  analysis  of  the  unloading  portion 
of  dynamic  HC  test  data  plotted  in  a  form  that  is 
conducive  to  constitutive  modeling. 
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Figure  3.48.  Results  from  the  COV  analysis  of  the  unloading  portions 
of  "measured"  dynamic  HC  test  data. 
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Figure  3.50.  Load-unload  responses  constructed  from  the  results  of 
the  COV  analysis  of  dynamic  HC  test  data  plotted  in  a 
form  that  is  conducive  to  constitutive  modeling. 
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Figure  3. 94.  Rosenblueth  equations  Tor  converting  the  results  from  the 
COV  analysis  of  "measured"  KQ  test  data  into  forms  that 
are  conducive  to  constitutive  modeling. 
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Figure  3.95.  Rosenblueth  equations  for  converting  the  results  from  the 
COV  analysis  of  "measured"  EXT  Kc  test  data  into  forms 
that  are  conducive  to  constitutive  modeling. 
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Figure  3.96.  Rosenblueth  equations  for  converting  the  results  from 
the  COV  analysis  of  "measured"  IIC  test  data  into  a 
form  that  is  conducive  to  constitutive  modeling. 


138 


PSD  = 


L/Ao 


S, 


S, 


e 

z 


PSD 


■(' 


2 


L/Ao 


Two  Variables 

for  Given  L/Ao  :  e  and  e 

z  r 


where  e  is  the 
rT 

true  radial  strain 


Two  Variables 

for  Given  L/Ao  :  e  and  e 

z  r 


e  is  er  *  (ep  ,  cr  +  oe  ,  er  -  ae  ) 

T  r*  p  ' 


where 


E  e 

—  r  start  of  shear 

er  N 

E(er  start  of  shear  er 

°e  N  -  1 

r 


ECS]  =  {[(S++)(1+p)  +  (S+")(1-p)  +  (s"+)(1-p)  +  (S  )(1+p)] 

V[S]  =  ^-[(S++)2(-p)  +  (S+")2(1-p)  +  (S_+)2(1-p)  +  (S'")2(1-p)] 

-  (E[S])2 

Figure  3.97.  Rosenblueth  equations  for  converting  the  results  from  the 
COV  analysis  of  "measured"  TXC  test  data  into  forms  that 
are  conducive  to  constitutive  modeling. 
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Figure  3.109.  Results  from  applying  the  Rosenblueth  procedure  to  the  results  from  the  COV  analysis  of 
"measured"  TXC  test  data  for  o  =0  and  0.1  MPa. 
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Figure  3.118.  Effect  of  including  p  in  the  Rosenblueth  procedure  for  TXC  test  data  at  o  =  17.5  MPa 
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Figure  3.122.  Results  from  applying  the  Rosenblueth  procedure  to 

constructed  continuous  results  from  the  COV  analysis  of 
"measured"  static  HC  test  data. 
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Figure  3.123.  Results  from  applying  the  Rosenblueth  procedure  to 

constructed  continous  results  from  the  COV  analysis  of 
"measured"  dynamic  11C  test  data. 


dU  'SS3H1S 

ItflXtf 

£> 

*dU  ’3DNja3JJlO  SS3«1S  lUdlDMiad 

1 — 1 
0 

QJ 

w 

3 

<b 

<D 

w 

C3 

\ 

0) 

2  2 

\ 

d) 

_B 

CE  C 

u  u 

\ 

u 

** 

2 

:  r 

\  \ 

£= 

U-J 

o 

1 

\  \ 

O 

1 

cn 

3 

\  ' 

U-< 

(/) 

•H 

w 

\  \  \ 

U) 

1 

u 

XJ 

<D 

• 

<u 

r-H 

1 

> 

\  \  \  \ 

CO 

u 

<3 

o  fft 

x  V  \  " 

□  § 

o 

c 

SS3ULS  idixu 


»d«  '33N3tf3Hl3  SSDblS  IddlDMlHd 


IS. 


cted  f 
V  anal 


MEAN  NORMAL  STRESS,  MPa 


-  MERN  ±  ? 


VOLUMETRIC  STRAIN  (cylinder),  PERCENT 


Figure  3.127.  Load-unload-reload  respi  uses  constructed  from  the 

results  of  applying  the  Rosenblueth  procedure  to  the 
results  from  the  COV  analysis  of  "measured"  static  HC 
test  data. 
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Figure  3.128.  Load-unload  responses  constructed  from  the  results  of 
applying  the  Rosenblueth  procedure  to  the  results  from 
the  COV  analysis  of  "measured"  dynamic  HC  test  data. 
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Figure  3. 129. 
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Figure  3. 131 


Comparison  of  the  results  from  the  COV  and  COV/ 
Rosenblueth  analysis  techniques  for  the  loading  portion 
of  static  HC  test  data. 
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Figure  3.132.  Comparison  of  the  results  from  the  COV  and  COV/ 

Rosenblueth  analysis  techniques  for  the  loading  portion 
of  dynamic  KC  test  data. 
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Figure  3  135.  Comparison  of  the  results  from  the  COV  and  COV/Rosenblueth  analysis  techniques  for  the 
loading  portion  of  static  TXC  test  data  at  o  =  0.7,  1.7,  and  3.5  MPa. 
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Figure  3.136.  Comparison  of  the  results  from  the  COV  and  COV/Rosenblueth  analysis  techniques  for 
the  loading  portion  of  static  TXC  test  data  at  a  =  17.5  MPa. 
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the  loading  portion  of  static  TXC  test  data  at  a  =  55.0  MPa. 
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Figure  3.140.  Comparison  of  the  results  from  the  COV  and  COV/Rosenblueth  analysis  techniques  for 
the  loading  portion  of  static  TXC  test  data  at  o  =  60.0  MPa. 


Figure  3.141.  Comparison  of  the  results  from  the  COV  and  COV/Rosenblueth  analysis  techniques  for 
the  loading  portion  of  static  TXC  test  data  at  a  =  80.0  MPa. 


CHAPTER  4 


FINAL  PROBABILISTIC  MECHANICAL  PROPERTIES 

4.1  APPLICABILITY  OF  THE  PROBABILISTIC  MECHANICAL  PROPERTIES  TO  AREA  A 

The  applicability  of  the  probabilistic  mechanical  properties  from 
Chapter  3  to  the  first  dry  density-based  layer  of  the  subsurface  profile 
for  Area  A  is  assessed  by  conducting  one-way  analysis  of  variance  tests 
on  the  dry  density  data  for  the  test  specimens  versus  the  dry  density 
data  for  the  layer.  The  results  of  the  analysis  for  each  test  type  are 
shown  in  Table  4.1.  The  mean  dry  density  for  the  test  specimens  used  in 
the  dynamic  uniaxial  strain  (UX)  and  hydrostatic  compression  (HC)  tests 
and  static  UX  with  lateral  stress  measurements  (Ko),  HC,  and  triaxial 
compression  (TXC)  tests  are  not  significantly  different  from  the  mean 
dry  density  for  the  layer.  The  mean  dry  density  for  test  specimens  used 
in  the  Ko  tests  conducted  in  the  extension  test  device  (EXT  Ko)  is 
statistically  different  from  the  mean  dry  density  for  the  layer.  There¬ 
fore,  except  for  the  EXT  Ko  test  data,  probabilistic  mechanical  proper¬ 
ties  based  directly  on  the  results  from  the  analysis  of  the  test  data 
are  assumed  to  be  adequate  for  the  layer.  The  results  from  the  analysis 
of  the  EXT  Ko  test  data  will  be  used  only  as  a  guide  for  developing  the 
negative  principal  stress  difference  (PSD)  regime  of  the  Ko  3tress  path 
unloading. 

4.2  PROBABILISTIC  MECHANICAL  PROPERTIES  DEVELOPED  FROM  TEST  RESULTS 

In  Section  3.2.5,  it  was  shown  that  the  results  from  applying 
covariance  analysis  to  test  data  plotted  in  forms  that  are  conducive  to 
constitutive  modeling  (COV  technique)  and  the  results  from  applying  the 
Rosenbluetn  procedure  to  the  results  from  the  covariance  analysis  of 
"measured  data"  (COV/Rosenblueth  technique)  were  essentially  the  same. 
Hence,  the  recommended  probabilistic  laboratory  responses  will  be  based 
on  the  results  from  the  COV  analysis  technique.  Unrealistic  responses 
from  the  analysis  of  the  laboratory  test  data  are  not  included  in  the 
recommended  responses. 

The  recommended  probabilistic  laboratory  UX  compressibility  and  the 
results  from  the  COV  analysis  technique  are  plotted  in  Figure  4.1  as 
thick  and  thin  lines,  respectively.  Since  the  axes  form  limits  to  the 
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possible  loading  UX  compressibility  response,  the  mean  ±  1  standard 

A  /\ 

deviation  (o)  responses  are  constructed  such  that  the  mean  -  3o  values 
are  0  or  positive.  The  recommended  probabilistic  response  does  not 
include  the  unrealistic  responses  from  the  COV  analysis,  i.e.,  the 
softening  of  the  loading  response  above  an  axial  stress  of  approximately 
80  MPa  and  the  unusual  "hook"  in  the  unloading  response  that  begins  at 
an  axial  stress  of  approximately  95  MPa. 

In  Figure  M.2,  the  recommended  probabilistic  laboratory  loading  UX 
stress  path  (thick  line)  is  drawn  slightly  below  the  results  from  the 
analysis  of  the  Ko  test  data  (thin  line).  From  Figure  3.15,  some  of  the 
test  data  shifted  in  an  increasing  principal  stress  difference  (PSD) 
direction  at  values  for  mean  normal  stress  (MNS)  of  approximately  25  and 
55  MPa.  The  unloading  and  reloading  stress  path  responses  are  con¬ 
structed  by  drawing  the  mean  response  first.  Then,  the  mean  ±  a  unload¬ 
ing  responses  are  drawn  by  laterally  shifting  the  mean  unloading  to 

A 

begin  at  the  appropriate  point  on  the  mean  ±  a  loading  response.  The 
response  is  extended  or  shortened  as  necessary.  From  Section  3. 2. 3. 2, 
there  is  very  little  variance  in  the  unloading  or  reloading  Ko  stress 
paths  from  a  given  stress  level.  The  reloading  responses  are  con¬ 
structed  in  the  same  way  but  begin  at  the  appropriate  point  on  the 
corresponding  unloading  response.  The  results  from  the  probabilistic 
analysis  of  the  second  set  of  EXT  Ko  unloading  responses  is  used  as  a 
guide  for  developing  the  negative  PSD  region  of  the  recommended  response 
(Figure  JJ.3). 

The  recommended  probabilistic  laboratory  static  HC  compressibility 
response  is  slightly  softer  and  has  a  greater  variance  than  the  probabi¬ 
listic  response  from  the  COV  analysis  of  the  static  HC  test  data 
(Figure  4.1J).  The  recommended  response  is  the  same  as  the  results  from 
the  COV  analysis  of  the  dynamic  HC  test  data  (Figure  ^.5).  Since  the 
results  from  the  COV  analysis  of  static  HC  test  data  are  within  the 

A 

mean  ±  o  response  for  the  dynamic  HC  test  data,  there  is  no  apparent 
dynamic  effect.  This  material  does  exhibit  a  slight  rate  effect  in  UX 
compressibility  for  static  to  dynamic  times  to  peak  stress  of  several 
minutes  to  less  than  10  milliseconds.  Therefore,  the  recommended  proba¬ 
bilistic  laboratory  HC  response  is  applicable  to  static  loadings  but  may 
not  be  applicable  to  fast  dynamic  loadings. 
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The  results  from  applying  the  COV  analysis  technique  to  the  loading 
static  TXC  shear  responses  at  constant  confining  pressures  less  than 
10  MPa  are  used  as  tne  recommended  probabilistic  laboratory  responses 
(Figures  4.6  and  4.7).  The  recommended  responses  and  the  results  from 
the  COV  analysis  on  comparable  test  data  for  constant  confining  pres¬ 
sures  of  20.0,  40.0,  60.0,  and  80.0  MPa  are  shown  in  Figures  4.8-4.11, 
respectively.  Since  there  is  very  little  variation  in  the  unloading  and 
reloading  responses,  only  a  mean  recommended  unloading  that  begins  at 
the  end  of  the  mean  loading  response  is  shown  for  each  confining  pres¬ 
sure. 

Recommended  probabilistic  laboratory  "peak"  strength  envelopes  are 
presented  in  Figure  4.12.  The  envelopes  are  eyeball  fits  through  the 
"peak"  strength  data  from  Table  3.12.  Slight  undulations  in  the  enve¬ 
lopes  are  generated  by  the  spline  fits  through  the  digitized  curves. 

The  "peak"  strengths  from  the  probabilistic  shear  responses  presented  in 
Figures  4.6-4.11  are  plotted  as  solid  circles.  The  probabilistic  shear 
responses  for  some  confining  pressures  have  "peak"  strengths  that  do  not 
follow  a  mean  ±  no  criteria.  However,  all  of  the  "peak"  strengths  from 
the  shear  responses  are  close  to  the  probabilistic  strength  envelope 
except  for  the  "peak"  strengths  at  a  confining  pressure  of  80  MPa. 

The  probabilistic  shear  response  for  a  confining  pressure  of 
80.0  MPa  will  be  adjusted  based  on  the  probabilistic  "peak"  strength 
envelope.  The  mean  shear  response  is  developed  using  the  mean  value  for 
"peak"  strength  from  Figure  4.12b  as  a  target  value,  and  the  recommended 
responses  for  confining  pressures  of  20.0,  40.0,  and  60.0  MPa  as  a  guide 
to  the  shape  of  the  response.  The  variance  is  based  on  the  coefficient 
of  variation  for  the  tests  at  confining  pressures  of  17.5,  20.0,  and 
55.0  MPa  and  the  recommended  responses  at  20.0  and  60.0  MPa 
(Figures  4.13  and  4.14).  The  coefficient  of  variation  starts  at  a  value 
of  approximately  50  percent  and  decreases  to  a  fairly  constant  value  of 
10  percent.  Using  these  criteria,  a  new  recommended  shear  response  for 
a  confining  pressure  of  80.0  MPa  is  shown  in  Figure  4.15. 

An  estimated  in  situ  effect  will  be  added  to  the  recommended  labo¬ 
ratory  dynamic  UX  compressibility  response  using  an  undocumented  proce¬ 
dure  based  on  the  analyses  of  two  dynamic  in  situ  compressibility 
tests.  A  relationship  was  observed  between  axial  stress  and  the  ratio 


of  seismic  strain  to  in  situ  strain  (see  Figure  4.16).  The  seismic 
strain  is  calculated  as 


e 


s 


(1+u)  C 


(4.1) 


where  e  is  the  seismic  axial  strain,  a  is  the  axial  stress,  g  is 
acceleration  due  to  gravity,  is  the  dry  density,  w  is  the  water 
content,  and  C  is  the  seismic  compression  wave  velocity.  The  in  situ 
strain  is  calculated  as 


eIN  SITU 


^ STRAIN 


(*».2) 


where  eIK  SITU  is  the  in  situ  axial  strain  and  &STRAjf<  io  obtained 
from  Figure  4.16.  The  Rosenblueth  procedure  is  used  since  probabilistic- 
values  for  dry  density  (1i4.7  ±  5.1  pcf)  and  water  content  (3.8  ± 

1.6  percent)  were  developed  in  Section  3.1.  Mean  and  standard  deviation 
values  for  C  of  2675  and  100  fps,  respectively,  are  typical  for  the 
area.  The  probabilistic  in  situ  dynamic  UX  compressibility  response  is 
shown  as  a  thick  line  in  Figure  4.3.  Since  the  procedure  for  calculat¬ 
ing  the  in  situ  axial  strain  is  only  valid  to  an  axial  stress  of  8  MPa, 
the  response  above  8  MPa  is  connected  to  the  recommended  laboratory 
response  (thin  line)  using  a  French  curve  and  the  author's  judgment. 

Methods  for  determining  the  in  situ  and  static  to  dynamic  effects 
for  the  UX  stress  path  of  this  material  are  not  available  at  this 
time.  Small  variations  in  the  UX  compressibility  response  do  not  appear 
to  significantly  effect  the  UX  stress  path.  Therefore,  the  in  situ  and 
dynamic  effects  incorporated  into  the  UX  compressibility  response  will 
probably  have  little  effect  on  the  UX  stress  path.  There  was  no  signif¬ 
icant  difference  observed  between  the  results  from  the  available  static 
and  dynamic  TXC  tests. 

4  3  METHODS  FOR  EXTRAPOLATING  THE  VARIANCE  OF  MECHANICAL  PROPERTIES  TO 

HIGHER  STRESSES 

Ground  shock  calculations  require  mechanical  properties  of  the 
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earth  materials  to  stresses  higher  than  those  generally  applied  in  the 
laboratory  tests.  Two  methods  for  extrapolating  the  variance  of  the 
recommended  responses  to  higher  stresses  will  be  examined.  The  mean 
nigh  pressure  responses  were  extrapolated  using  engineering  knowledge. 

The  coefficient  of  variation  and  values  of  alpha  (a)  and  beta  (B) 
for  a  beta  distribution  will  be  calculated  for  the  recommended  probabi¬ 
listic  laboratory  UX  compressibility,  UX  stress  path,  and  TXC  "peak” 
strength  envelope.  The  trends  for  these  parameters  will  be  used  to 
determine  the  variance  about  the  extrapolated  mean  responses.  The 
coefficient  of  variation  is  calculated  as  the  standard  deviation  divided 
by  the  mean  and  is  presented  as  a  percentage.  The  a  and  B  parame¬ 
ters  for  the  beta  distribution  are  calculated  as 

2 

a  =  —  (1  -  x)  -  (1  +  x)  (H.3) 

V 


B  =  -  (a  +  2)  (HJD 

x 

—  x  ~  a  — 

where,  x  =  - -  and  V  = 

b  -  a 

These  equations  are  based  on  the  assumption  that  the  absolute  minimum 
(  a  )  and  maximum  (  b  )  values  are  zero  and  mean  *  3 a  ,  respectively. 

Plots  of  coefficient  of  variation,  a  ,  and  B  versus  stress  are 
shown  in  Figures  4. 18-4.23.  The  values  for  the  recommended  probabilis¬ 
tic  laboratory  responses  and  the  results  from  the  COV  analysis  of  the 
test  data  are  shown  as  thick  and  thin  lines,  respectively.  The  plots 
tend  to  become  parallel  to  the  X-axis  as  stress  increases.  Values  for 
coefficient  of  variation,  a  ,  and  B  of  20  percent,  8.5,  and  5,  respec¬ 
tively,  for  the  UX  compressibility,  10  percent,  22,  and  6.7,  respective¬ 
ly,  for  the  UX  stress  path,  and  5  percent,  55,  and  7.8,  respectively  for 
the  TXC  "peak"  strength  envelope  will  be  used. 

The  extrapolated  responses  are  shown  in  Figures  4. 24-4. 26.  For  the 
UX  compressibility  and  stress  path,  extrapolating  the  variance  using 
coefficient  of  variation  and  the  beta  distribution  parameters  provide 
the  same  mean  ±  o  responses;  therefore,  coefficients  of  variation  of  20 
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and  10  percent  will  be  used  to  extrapolate  the  variance  to  higher 
stresses  since  it  is  the  simplest  to  apply  and  requires  no 
assumptions.  For  the  TXC  "peak"  strength  envelope  (Figure  4.26),  a 
coefficient  of  variation  of  5  percent  will  be  used  since  it  provided  a 
slightly  larger  variance  than  the  beta  distribution  parameters. 

4.4  FINAL  PROBABILISTIC  MATERIAL  PROPERTIES 

The  low  stress  portion  of  the  recommended  probabilistic  dynamic  UX 
compressibility  and  stress  path,  laboratory  HC  compressibility,  TXC 
"peak"  strength  envelope,  and  TXC  shear  responses  are  shown  in  Figures 
4.27-4.32.  They  are  the  same  as  those  shown  in  Figures  4.1-4.10,  4.12, 
4.15,  and  4.17.  A  deterministic  triaxial  extension  "peak"  strength 
envelope  is  shown  in  Figure  4.29.  The  HC  response  will  not  be  extrapo¬ 
lated  to  higher  stresses  since  there  is  no  available  data  to  use  as  a 
guide. 

A  major  concern  with  developing  the  probabilistic  high  stress 
properties  is  the  correct  handling  of  the  full  saturation  portion  of  the 
responses,  i.e.,  when  the  air  voids  content  is  zero  and  the  degree  of 
saturation  is  100  percent.  The  UX  compressibility  response,  UX  stress 
path,  and  TXC  "peak"  strength  envelope  will  have  slopes  that  are  nearly 
constant  for  stress  greater  than  or  equal  to  the  stress  at  full  satura¬ 
tion  (the  "lockup"  stress).  First  consider  the  results  of  applying  the 
covariance  analysis  technique  to  a  set  of  test  data  that  achieve  full 
saturation  at  different  stress  levels.  Since  the  slope  of  the  responses 
is  analyzed,  a  constant  slope  for  the  probabilistic  response  will  not  be 
reached  until  all  of  the  responses  in  the  data  set  have  a  constant 
slope.  This  indicates  that  the  "lockup"  stress  for  the  probabilistic 
response  is  the  highest  "lockup"  stress  for  the  data  set  and  is  not 
probabilistic. 

The  high  stress  probabilistic  dynamic  in  situ  UX  compressibility 
response  is  shown  in  Figure  4.33*  Below  an  axial  stress  of  93  MPa,  the 
response  is  the  same  as  that  shown  in  Figure  4.27.  The  mean  ±  0  re¬ 
sponse  above  an  axial  stress  of  93  MPa  was  determined  by  multiplying  the 
mean  axial  strain  at  a  given  value  of  axial  stress  by  a  coefficient  of 
variation  of  20  percent.  The  mean  -  j  ,  mean,  and  mean  +  0  responses 
have  constant  moduli  of  50,  45,  and  40  GPa  above  a  "lockup"  axial  stress 
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of  700  MPa.  The  mean  unloading  response  from  an  axial  stress  of  700  MPa 
was  determined  using  engineering  judgment.  From  the  lower  pressure 
probabilistic  unloading  responses,  the  coefficient  of  variation  for 
axial  strain  at  the  end  of  an  unloading  set  is  20  to  23  percent.  A 
coefficient  of  variation  of  23  percent  is  used  to  determine  the  strain 
at  zero  axial  stress  for  the  final  mean  ±  o  unloadings. 

The  high  stress  probabilistic  dynamic  in  situ  UX  stress  path  and 
TXC  "peak"  strength  envelope  are  shown  in  Figure  4.3*1.  A  deterministic 
triaxial  extension  "peak"  strength  envelope  is  also  shown.  The  UX 
stress  path  below  a  MNS  of  70  MPa  is  the  same  as  that  shown  in 
Figure  4.29;  the  TXC  "peak"  strength  envelope  below  a  MNS  of  1 60  MPa  is 
the  same  as  that  shown  in  Figure  4.14.  Values  for  coefficient  of  varia¬ 
tion  of  10  and  5  percent  were  used  to  extrapolate  the  variance  of  the  UX 
stress  path  and  TXC  "peak"  strength  envelope,  respectively.  Constant 
loading  and  unloading  slopes  that  imply  values  for  Poisson's  ratio  of 
0.487,  0.490,  and  0.493  are  used  for  the  mean  -  o  ,  mean,  and 
mean  +  a  UX  stress  paths,  respectively,  for  MNS  greater  than  or  equal  to 
572  MPa.  The  probabilistic  TXC  "peak"  strength  envelope  becomes  flat  at 
a  MNS  of  500  MPa.  On  the  mean  responses,  these  values  of  MNS  coincide 
with  the  "lockup"  axial  stress  of  700  MPa.  The  mean  UX  stress  path 
unloadings  are  developed  using  engineering  judgment.  A  coefficient  of 
variation  of  10  percent  was  used  to  determine  the  initial  slopes  for  the 
mean  ±  o  unloadings.  The  unloadings  were  completed  by  attaching  the 
mean  unloading  response  to  the  initial  portions. 
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Figure  4.1.  Comparison  of  the  results  from  the  COV  analysis  of  dynamic 
UX  test  data  and  the  recommended  laboratory  response . 
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Figure  4.4.  Comparison  of  the  re; 

HC  test  data  and  the 
bility  response. 


ults  from  the  COV  analysis  of  static 
recommended  laboratory  HC  compressi- 


MERN  NORMRL  STRESS 


Figure  4.6.  Results  from  the  COV  analysis  of  the  static  TXC  test  data  and  the  recommended  unloading 
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Figure  4.9.  Comparison  of  cho  results  from  the  COV  analysis  of  the  static  TXC  test  data  and  the  recommended 
laboratory  response  for  a _  =  40.0  MPa. 
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Figure  4.13.  Coefficient  of  variation  for  the  results  from  the  COV  analysis  of  static  TXC  test  data 
for  o  =  17.5  and  20.0  MPa  and  the  recommended  laboratory  response  for  o  =  20.0  MPa 
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Figure  4.16.  Correlation  between  in  situ  axial  strain  and  seismic 
axial  strain. 
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Figure  4.17.  Recommended  laboratory  and  in  situ  dynamic  UX 
compressibility  responses. 
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Figure  A. 18. 


Coefficient  of  variation  for  the  results  from  the  COV 
analysis  of  dynamic  UX  test  data  and  the  recommended 
laboratory  UX  compressibility  response. 


Figure  4.20.  Coefficient  of  variation  for  the  results  from  the  COV  analysis  of  K  and  EXT  K 
data  and  the  recommended  laboratory  stress  path.  °  ' 
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Figure  4.21.  Beta  distribution  parameters  a  and  B  lor  the  results 
from  the  COV  analysis  of  KQ  and  EXT  K0  test  data  and  the 
recommended  laboratory  stress  path. 
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Figure  A. 27.  Low-stress  laboratory  and  in  situ  probabilistic  dynamic  UX 
compressibility  responses. 
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strength  envelope  and  probabilistic  TXC  "peak"  strength 
envelope  and  UX  stress  path. 


CHAPTER  5 


SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 

5.1  SUMMARY 

A  procedure  for  developing  probabilistic  properties  of  earth  mate¬ 
rials  for  use  in  probabilistic  ground  shock  calculations  was  presented 
in  Chapter  2.  In  this  procedure,  the  subsurface  profile  of  a  given  area 
was  based  on  material  type,  mass- volume  relations,  index  values,  and 
seismic  velocity.  Probabilistic  values  for  the  mass-volume  relations, 
index  parameters,  and  seismic  velocity  were  determined  using  simple 
statistical  procedures.  Probabilistic  mechanical  properties  (i.e., 
strength  and  compressibility  relations)  were  developed  based  on  the 
results  of  laboratory  tests  conducted  on  undisturbed  samples  from  each 
layer  of  the  subsurface  profile  and  on  field  experiments  and  engineering 
judgment.  Data  obtained  from  two  adjacent  areas  located  near  Yuma, 
Arizona  were  used  to  illustrate  the  application  of  the  procedure. 

Because  the  basic  material  type  at  the  two  areas  does  not  change 
significantly  in  the  upper  150  feet,  the  subsurface  profile  was  based  on 
mass-volume  relations,  index  parameters,  and  seismic  velocity.  Several 
probable  zones  were  identified  from  plots  of  dry  density  data  versus 
depth  (Section  3.1).  The  apparent  difference  in  the  data  was  evaluated 
by  conducting  one-way  analysis  of  variance  tests.  If  the  test  results 
indicated  no  statistically  significant  difference,  the  trends  were 
reevaluated  and  the  test  repeated.  Probabilistic  values  for  wet  den¬ 
sity,  degree  of  saturation,  and  air  voids  content  were  calculated  using 
a  procedure  developed  by  Rosenblueth  (Section  2.2.3). 

The  results  from  uniaxial  strain  (UX),  UX  with  lateral  stress 
measurement  (Ko) ,  hydrostatic  compression  (HC),  and  triaxial  compression 
(TXC)  tests  conducted  on  undisturbed  samples  from  the  first  dry  density- 
based  layer  at  Area  A  were  used  to  illustrate  the  procedure  for  devel¬ 
oping  probabilistic  mechanical  properties.  The  data  base  was  supple¬ 
mented  by  tests  on  undisturbed  samples  obtained  from  the  upper  ^9  feet 
of  Area  B  and  with  values  of  dry  density  within  one  standard  deviation 
of  the  mean  dry  density  for  Area  A. 
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A  procedure  for  removing  electrical  noise  from  the  data 
(Section  2.2.1)  was  applied  to  most  of  the  UX  test  data  and  the  higher 
pressure  Ko  and  TXC  test  data.  The  mean  and  variance  of  incremental 
slopes  from  plots  of  the  smoothed  data  were  analyzed  using  the  covari¬ 
ance  technique  (Section  2.2.2).  The  covariance  (COV)  analysis  was 
conducted  on  (a)  smoothed  data  plotted  in  a  way  that  was  conducive  to 
constitutive  modeling,  i.e.,  axial  stress,  mean  normal  stress,  principal 
stress  difference,  axial  strain,  volumetric  strain,  and  principal  strain 
difference,  and  (b)  smoothed  data  plotted  in  a  way  that  was  as  close  as 
practical  to  the  responses  actually  measured  during  the  test,  i.e., 
stress,  load,  and  deflection.  The  results  from  the  second  analysis  were 
converted  into  the  conducive  forms  for  consititutive  modeling  by  apply¬ 
ing  the  Rosenblueth  procedure.  The  load,  unload,  and  reload  portions  of 
the  data  were  analyzed  separately  since  the  COV  analysis  requires  that 
the  individual  responses  be  monotonic. 

In  the  Rosenblueth  procedure  (Section  3.2.JJ),  the  mean  and  mean  ±  1 
standard  deviation  (a)  secant  slopes  were  calculated  using  the  results 
from  the  COV  analysis  of  the  "measured"  mechanical  property  test  data. 
The  probabilistic  target  responses  were  constructed  by  multiplying  these 
values  by  the  mean  value  for  the  denominator  of  the  slope.  Since  the 
slopes  of  the  desired  plots  were  calculated,  it  was  necessary  to  inves¬ 
tigate  the  effect  of  including  correlation  coefficients  (p).  The  values 
for  p  were  estimated  by  linear  regression  analysis  conducted  on  the 
random  variables.  Including  p  had  the  greatest  effect  on  the  variance 
and  little  to  no  effect  on  means.  Whether  the  variance  was  increased, 
decreased,  or  unchanged  depended  on  the  sign  of  p  and  the  rank  of  the 
++,  -+,  and  —  calculations. 

The  results  from  applying  the  two  analysis  techniques  to  loading 
test  data  were  compared  in  Section  3-2.5.  Both  techniques  provided 
similar  results  for  the  Ko,  EXT  Ko,  and  HC  test  data.  For  TXC  test 
data,  the  results  were  slightly  different  since  the  stress-axis  was 
incremented  in  the  COV  analysis  of  the  "measured  data."  The  unloading 
and  reloading  responses  from  both  analysis  techniques  were  essentially 
the  same  for  all  test  types. 

Low  stress  probabilistic  mechanical  properties  for  the  layer  were 
developed  by  using  the  results  from  the  COV  analysis,  analyses  of 


dynamic  in  situ  compressibility  tests,  and  engineering  judgment 
(Section  4.2).  The  variance  of  the  low  stress  probabilistic  UX  com¬ 
pressibility,  UX  stress  path,  and  TXC  "peak"  strength  envelope  was 
extrapolated  to  higher  stresses  by  examining  the  coefficients  of  varia¬ 
tion  and  beta  distribution  parameters  for  the  lower  stress  responses 
(Section  4.3);  constant  values  were  approached  as  stress  increased. 

The  final  probabilistic  mechanical  properties  for  the  layer  were 
presented  in  Section  4.4.  The  low  stress  responses  were  extended  past 
the  stress  level  at  which  the  material  became  fully  saturated  ("lockup" 
stress).  The  "lockup"  stress  was  not  treated  probabilistically  since  a 
single  value  would  have  resulted  from  the  COV  analysis  of  a  set  of  test 
results  that  achieved  different  values  for  "lockup"  stress. 

5.2  CONCLUSIONS 

The  following  conclusions  can  be  made  as  a  result  of  this 
investigation: 

1.  The  one-way  analysis  of  variance  test  is  very  useful  for  devel¬ 
oping  the  subsurface  profile  of  a  given  area  based  on  plots  of  mass- 
volume,  index,  and  seismic  velocity  data  versus  depth. 

2.  The  COV  analysis  of  mechanical  property  test  data  for  earth 
materials  plotted  in  fo'mns  that  are  conducive  to  constitutive  modeling 
and  applying  the  Rosenblueth  procedure  to  the  results  from  the  COV 
analysis  of  "measured  data"  provide  similar  probabilistic  responses. 

3.  The  unioueness  of  a  soil  specimen's  response  during  a  mechan¬ 
ical  property  test  must  be  maintained  when  plotting  measured  data.  This 
is  accomplished  by  normalizing  the  measured  data  based  on  the  specimen 
dimensions.  A  hypothetical  example  illustrating  the  effect  of  using  the 
measured  deflection  is  shown  in  Figure  5.1.  Assuming  that  the  variance 
in  the  strain  of  two  test  specimens  of  different  size  is  the  result  of 
random  scatter  and  not  size  effect,  the  average  strain  should  be 

12.5  percent  and  not  13.0  percent  as  calculated  using  the  measured 
del  lection. 

4.  When  using  the  Rosenblueth  procedure  to  transform  the  results 
from  the  COV  analysis  of  the  "measured  data"  into  the  forms  that  are 
more  conducive  to  constitutive  modeling,  the  effect  of  correlation 
coefficients  p  must  be  evaluated. 


5.  The  procedure  used  fco  extrapolate  the  laboratory-based 
probabilistic  mechanical  properties  maximized  the  variance  about  the 
extrapolated  mean.  Tests  to  higher  stresses  are  needed  to  substantiate 
the  extrapolation  procedure. 

6.  The  mean  values  for  the  mass-volume  relations,  index  parame¬ 
ters,  UX  compressibility,  UX  stress  path,  TXC  shear  responses,  and  TXC 
"peak"  strength  envelope  are  coupled,  i.e.,  they  would  be  used  in  a 
deterministic  ground  shock  calculation  to  obtain  mean  values  for  the 
output.  Since  there  was  no  apparent  tie  between  the  mean  ±  a  proper¬ 
ties,  the  mean  *  c  or  mean  -  a  values  for  one  material  property  cannot 
be  coupled  with  the  mean  +  o  or  mean  -  a  value  for  the  other  material 
properties.  This  results  in  many  possible  combinations. 

5.3  RECOMMENDATIONS 

Recommendations  for  further  study  are: 

1.  Include  SDatial  variability  in  the  determination  of  probabi¬ 
listic  properties  of  earth  materials. 

2.  Correlate  the  probabilistic  mass-volume,  index,  and  mechanical 
properties  of  earth  materials  with  the  results  from  easily  obtained 
field  tests  such  as  cone  penetration  tests. 

3.  Examine  the  true  distribution  of  mass-volume,  index,  and  me¬ 
chanical  property  data  to  quantify  the  effect  of  assuming  a  symmetrical 
distribution  in  the  analysis  technique. 

4.  Complete  the  analysis  for  all  layers  of  the  profile  for  a  given 
site  and  determine  if  the  differences  in  the  mass-volume,  index,  and 
seismic  velocity  data  is  translated  into  different  mechanical 
responses.  A  method  such  as  one-way  analysis  of  variance  could  be  used 
to  determine  if  the  apparent  differences  are  significant. 


Sample  A  is  2"  high 
Sample  B  is  3"  high 

Assume  that  during  a  UX  strain  test, 

at  o  =20  MPa  , 
z 

z.  =  102  and  e  =  152 
a  B 

The  average  strain  is  1 — —  =  12.52 

AHa  =  2(0.10)  =  0.2" 

AHd  =  3(0.15)  =  0.45" 

D 

02+0  45 

The  average  AH  =  — 1 — - — : —  =  0.325" 
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The  average  H  =  — - —  =  2.5" 

The  average  strain  using  the  average  AH 
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Figure  5.1.  Discrepancy  created  by  varying  sample  size. 


REFERENCES 


1.  Chou,  Yu  T«  September  1981.  "Probabilistic  Modeling  of  Sub¬ 
grade  Soil  Strengths,"  Miscellaneous  Paper  GL-81-5,  US  Army  Engineer 
Waterways  Experiment  Station,  Vicksburg,  MS. 

2.  Diaz-Padiila,  Jorge,  and  Vanmarcke,  Erik  H.  August  1975. 
"Settlement  Prediction:  A  Probabilistic  Approach,"  Research  Report  R73~ 
40,  Department  of  Civil  Engineering,  Massachusetts  Institute  of  Technol¬ 
ogy,  Cambridge,  MA. 

3.  Gilbert,  Lawrence  W.  July  1977.  "A  Probabilistic  Analysis  of 
Embankment  Stability  Problems,"  Miscellaneous  Paper  S— 77— 1 0,  US  Army 
Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 

4.  Vanmarcke,  Erik  H.  "Reliability  of  Earth  Slopes,"  Journal  of 
the  Geotechnical  Engineering  Division,  Proc.  ASCE,  Vol.  103,  No.  GT1 1 , 
November  1977,  pp  1247-1265. 

5.  Jackson,  J.  G.  July  1982.  "Site  Characterization  for  Probabi¬ 
listic  Ground  Shock  Predictions,"  Miscellaneous  Paper  SL-82-8,  US  Army 
Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 

6.  Vanmarcke,  Erik  H.  December  1979.  "Probabilistic  Soil  Sam¬ 
pling  Program  for  MX-Related  Site  Characterization,"  Miscellaneous  Paper 
SL-79-26,  US  Army  Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 

7.  Vanmarcke,  Erik  H.  "Probabilistic  Modeling  of  Soil  Profiles," 
Journal  of  the  Geotechnical  Engineering  Division,  Proc.  ASCE,  Vol.  103, 
No.  GT1 1 ,  November  1977,  pp  1227-1246. 

8.  Benjamin,  Jack  R.,  and  Cornell,  C.  Allen.  "Probability,  Sta¬ 
tistics,  and  Decision  for  Civil  Engineers,"  McGraw-Hill  Book  Company, 
1970. 

9.  Rohani,  Behzad.  June  1982.  "Probabilistic  Solution  for  One- 
Dimensional  Plane  Wave  Propagation  in  Homogeneous  Bilinear  Hysteretic 
Materials,"  Miscellaneous  Paper  SL-82-7,  US  Army  Engineer  Waterways 
Experiment  Station,  Vicksburg,  MS. 

10.  Rohani,  Behzad,  and  Cargile,  James  D.  April  1984.  "A  Probabi¬ 
listic  One-Dimensional  Ground-Shock  Code  for  Layered  Nonlinear  Hyste¬ 
retic  Materials,"  Miscellaneous  Paper  SL-84-6,  US  Army  Engineer  Water¬ 
ways  Experiment  Station,  Vicksburg,  MS. 

11.  Baladi,  George  Y.,  and  Barnes,  Donald  E.  October  1983.  "An 
Objective  Waveform  Comparison  Technique,"  Technical  Report  SL-83-4,  US 
Army  Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 


12.  Vanmarcke,  Erik  H.,  Jackson,  Andrew  E.  Jr.,  and  Akers,  Stephen 
A.  September  1986.  "Ralston  Valley  Soil  Compressibility  Study:  Corre¬ 
lation  of  Laboratory  Uniaxial  Strain  Compressibility  with  Field  Seismic 
Refraction  Data,"  Technical  Report  SL-86-1 4,  US  Army  Engineer  Waterways 
Experiment  Station,  Vicksburg,  MS. 

13-  Jackson,  A.  E.  Jr.  December  1984.  "Statistical  Variation  of 
Laboratory  Soil  Properties,"  Miscellaneous  Paper  SL-84-1 4,  US  Army 
Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 

14.  Miller,  Irwin,  and  Freund,  John  E.  "Probability  and  Statistics 
for  Engineers,"  Prentice-Hall,  Englewood  Cliffs,  New  Jersey,  1977. 

15.  Baladi,  G.  Y.,  and  Rohani,  B.  June  1987.  "A  Method  for  the 
Statistical  Analysis  of  the  Stress-Strain  Properties  of  Earth  Materi¬ 
als,"  Proceedings  of  the  Thirty-Second  Conference  on  the  Design  of 
Experiments,  Report  No.  87-2,  U.  S.  Army  Research  Office,  Research 
Triangle  Park,  North  Carolina. 

16.  Rosenblueth,  E.  October  1975.  "Point  Estimates  for  Probabil¬ 
ity  Moments,"  Proceedings  of  the  National  Academy  of  Sciences,  Vol  72, 
No.  10,  pp  3812-381  4. 

17.  Rosenblueth,  E.  October  1981.  "Two-Point  Estimates  in  Proba¬ 
bilities,"  Applied  Mathematical  Modeling,  Vol  5,  pp  329-335. 

18.  Bowles,  Joseph  E.  "Physical  and  Geotechnical  Properties  of 
Soils,"  McGraw-Hill  Book  Company,  1979. 

19.  Akers,  S.  A.,  Reed,  P.  A.,  and  Ehrgott,  J.  Q.  August  1986. 

"WES  High-Pressure  Uniaxial  Strain  and  Triaxial  Shear  Test  Equipment," 
Miscellaneous  Paper  SL-86-1 1,  US  Army  Engineer  Waterways  Experiment 
Station,  Vicksburg,  MS. 

20.  Ehrgott,  J.  Q.  May  1971.  "Calculation  of  Stress  and  Strain 
from  Triaxial  Test  Data  on  Undrained  Soil  Specimens,"  Miscellaneous 
Paper  S-71-9,  US  Army  Engineer  Waterways  Experiment  Station.  Vicksburg, 
MS. 


232 


DISTRIBUTION  LIST 


DEPARTMENT  OF  DEFENSE 
Director 

Defense  Nuclear  Agency 
ATTN:  DFTD  (Dr.  D.  A.  Linger) 

DFSP  (Dr.  G.  W.  Ullrich) 

SPWE  (Mr.  C.  B.  McFarland) 

SPSD  (Dr.  K.  L.  Goerlng) 

SPSD  (MAJ  M.  A.  Reed) 

SPWE  (MAJ  Mike  Pelkey) 

SPWE  (Dr.  E.  J.  Rinehart) 
Technical  Library 
Washington,  DC  20305-1000 

Director 

Defense  Nuclear  Agency 
Nevada  Operations  Office 
ATTN:  TDNV  (Mr.  J.  W.  LaCornb) 

P.0.  Box  98518 

Las  Vegas,  NV  89193-8518 

Defense  Nuclear  Agency 
ATTN:  TDTT 
TDST 

Kirtland  AFB,  NM  87115-5000 
Director 

Defense  Advanced  Research  project  Agency 
ATTN;  Technical  Library 
11400  Wilson  Blvd. 

Arlington,  VA  22209 

Director 

Defense  Intelligence  Agency 
ATTN:  Technical  Library 
Washington,  DC  20301-6111 

Defense  Tecnnlcal  Information  Center 
ATTN:  TC 
Cameron  Station 
Alexandria,  VA  2231*1 


DEPARTMENT  OF  THE  ARM* 

Commander 

US  Array  Corps  of  Engineers 
ATTN:  CERD-L 

CERD-M  (Mr.  B.  0.  Benn) 

CEEC-ET  (Mr.  R.  L.  Wight) 

CEIM-SL 

Washington,  DC  20311-1000 
Division  Engineer 

US  Array  Engineer  Division,  Huntsville 
ATTN:  CEHND-SR 
P.0.  Box  1600 

Huntsville,  AL  35807-1301 
Director 

US  Array  Construction  Engineering  Research 
Laboratory 

ATTN:  Technical  Library 
P.0.  3ox  1005 
Chaapaign,  IL  61820-1305 


DEPARTMENT  OF  THE  ARM*  (CONTINUED) 

District  Engineer 
US  Army  Engineer  District,  Omaha 
ATTN:  CEMRO-ED-S  (Mr.  Bob  Kelley) 
CEMRO-ED-SH  (Mr.  BUI  Gaube) 

215  N.  17th  Street 
Omaha,  NE  68102-1978 

Commander/Director 
US  Army  Cold  Regions  Research  and 
Engineering  Laboratory 
ATTN:  Technical  Library 
72  Lyme  Road 
Hanover,  NH  03755-1290 

Commandant 

US  Army  Engineer  School 
ATTN:  ATZA-CD  (COL  Fred  Parker) 
Technical  Library 
Fort  Belvoir,  VA  22060-5281 

Commander 

Harry  Diamond  Laboratories 
Department  of  the  Army 
ATTN:  Technical  Library 
2800  Powder  Mill  Road 
Adelphl,  MD  20783-1197 

Director 

US  Army  Ballistic  Research  Laboratory 

ATTN:  Technical  Library 

Aberdeen  Proving  Ground,  KD  21005-5066 

Commander 

US  Army  Nuclear  and  Chemical  Agency 
ATTN:  Technical  Library 
7500  Backlick  Road,  Bldg.  2073 
Springfield,  VA  22150 


DEPARTMENT  OF  THE  NAVY 

Naval  Civil  Engineering  Laboratory 
ATTN:  Technical  Library 
Port  Huenerae,  CA  93013 

Naval  Facilities  Engineering  Command 
200  Stoval  Street 
ATTN:  Technical  Library 
Alexandria,  VA  22332 


DEPARTMENT  OF  THE  AIR  FORCE 

Air  Force  Institute  of  Technology 
Air  University 
ATTN:  Technical  Library 
Wrlght-Patterson  AFB,  OH  15133 

Air  Force  Office  of  Scientific  Research 
ATTN:  Technical  Library 
Bolling  AFB,  DC  20332 


DEPARTMENT  OF  THE  AIR  FORCE  (CONTINUED) 


Air  Force  Weapons  Laboratory  (AFSC) 

ATTN:  NTESG  (CPT  C.  W.  Felice) 

NTEDC  (Mr.  J.  N.  Thomas) 

Technical  Library 
Xirtland  AFB,  NM  87117-6008 

Air  Force  Engineering  and  Services  Center 
( AFSC) 

ATTN:  Technical  l  .b^ary 
Tyndall  AFB,  FL  32403 

Commander 

Ballistic  Missile  Office  (AFSC) 

ATTN:  HYEB  (LTC  D.  H.  Gage) 

ENSS  (LT  K.  R.  Hlrlinger) 

Technical  Library 
Norton  AFB,  CA  92*109-61168 


DEPARTMENT  OF  ENERGY 

Lawrence  Livermore  National  Laboratory 
ATTN:  Technical  Library 
P.0.  Box  808 
Livermore,  CA  9*1550 

Los  Alamos  National  Laboratory 
ATTN:  Technical  Library 
P.0.  Box  ’663 
Los  Alamos,  NM  875*15 

Sandla  National  Laboratories 
ATTN:  Technical  Library 
P.O.  Box  5800 
Albuquerque,  NM  87185 

Sandla  National  Laboratories 
ATTN:  Technical  Library 
Livermore,  CA  9*1550 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Mr.  J.  L.  Bratton 
Applied  Research  Associates,  Inc. 
*1300  San  Mateo  Blvd.,  NE,  Suite  A220 
Albuquerque,  NM  87110 


Dr.  Hon-Ylm  Ko 

Department  of  Civil,  Environmental,  and 
Architectural  Engineering 
University  of  Colorado  at  Boulder 
Boulder,  CO  80309 

Dr.  K.  C.  Vala.nl s 
Endochronlcs,  Inc. 

8605  Northwest  Lakecrest  Court 
Vancouver ,  VA  98665 


New  Mexico  Engineering  Research  Institute 
ATTN:  Technical  Library 
University  of  New  Mexico 
Box  25,  University  Station 
Albuquerque ,  NM  871 31 

Mr.  J.  G>  Lewis 

Dr.  Don  Simons 

RAD  Associates 

P.O.  Box  9695 

Marina  del  Rey,  CA  90291 

Mr.  L.  S.  Melzer 

Science  Applications  International 
Corporation 
505  We3t  Texas  Street 
First  City  Center,  Tower  2,  Suite  1335 
Midland,  TX  79701 

Dr.  John  Schatz 

Science  Applications  International 
Corporation 
P.O.  Box  2351 
La  Jolla,  CA  92038-2351 

Dr.  H.  E.  Read 

Mr.  Steve  Peyton 

S-Cubed 

P.O.  Box  1620 

La  Jolla,  CA  92038-1620 

Dr.  Lytih  Seaman 
SRI  International 
333  Raver.swood  Avenue 
Menlo  Park,  CA  99025 

Mr.  S.  J.  Grsen 
Mr.  William  Klauber 
Terra  Tek,  Inc. 

420  Wakara  Way 

Salt  Lake  City,  UT  84108 


Mr.  S.  E.  Blouln 

Applied  Research  Associates,  Inc. 
Box  120A,  Waterman  Road 
South  Royalton,  VT  05068 

Mr.  J,  L.  Drake 

Applied  Research  Associates,  Inc. 
3202  Wisconsin  Avenue 
Vicksburg,  MS  39180 


Mr.  Norman  Lipner 

Dr.  M.  G.  Katona 

TRW  Defense  Systems  Group 

P.O.  Box  1310 

San  Bernardino,  CA  92402 


I 


Dr.  J.  0.  Trullo 
Applied  Theory,  Inc. 
930  S.  LaBrea  Avenue 
Los  Angeles,  CA  90036 


Dr.  7.  Marvin  Ito 
Mr.  S.  H.  Schuster 

California  Research  A  Technology,  Inc. 
20943  Devonshire  Street 
Chatsworth,  CA  91311-2376 


>  KV 


Dr.  D.  J.  Ness 

TRW  Defense  and  Space  Systems  Group 
One  Space  Park,  Bldg.  134/Rm  9835 
Redondo  Beach,  CA  90273 

Dr.  I.  S.  Sandler 
Weldlinger  Associates 
333  Seventh  Avenue 
New  York,  NY  10001 


234 


vw  v_  >r,  vv 


1  V  V\ 


Dr.  Jereay  Isenbsrg 
Dr.  Howard  S.  Levine 
Weidiinger  Associates 
620  Hansen  Hay,  Suite  100 
Palo  Alto,  C.1  9*130!; 


Dr.  H.  F.  Carroll 
Departaer.t  of  Civil  Engineering 
and  Environmental  Sciences 
University  of  Central  Florida 
Orlando,  FL  328! 6 


Dr.  R.  T.  Allen 
Pacifica  Technology 
P.0.  Box  1H5 
Del  Mar,  CA  920111 


Dr.  ?.  H.  S.  W.  Kuiati lake 
University  of  Arizona 
3uilding  ”o.  12 
Tuscan,  AZ  85721 


Dr.  Robe-'-.  M.  Schoites 
Dr.  Ralph  S-"no 
Mississippi  State  University 
°.0.  Draper  CE 

Mississippi  State,  MS  39762 


235 


t 

a 

V 


^  «na  *+.  -  i  «  >,  >,  «— ».«  M  *.  *  ..  a  ■«  m  ,,  *  *  %  «  *  m>*  n  *  m  --Q.w  ■*<*■■*  >.  A  ■*,>  —  A  *  .1  »  f,  .  f  ,fc  — 


